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ABSTRACT 


This  report  is  concerned  with  formulas  for  the  determination  of  the  earth’s 
physical  surface  and  external  gravity  field  from  free -air  gravity  anomalies  to 
an  approximation  linear  in  the  elevation  and  its  derivatives. 

Part  A  considers  integral  equations  and  their  linear  solutions;  Part  B  gives 
an  elementary  deduction  of  these  solutions  from  the  geometrically  evident 
gradient  solution;  and  the  subject  of  Part  C  is  an  application  to  various  gravity- 
dependent  quantities  and  an  evaluation  of  different  solutions. 
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INTRODUCTION 


The  present  work  is  a  continuation  of  an  earlier  report  (Moritz,  1964). 

The  previous  report  gave  a  discussion  of  the  various  formulations  of  the  geodetic 
boundary -value  problem,  that  is,  the  gravimetric  determination  of  the  geoid  or 
of  the  physical  surface  of  the  earth,  using  either  unreduced  free-air  anomalies 
or  other  anomalies  that  correspond  to  various  gravity  reductions.  The  present 
report  is  exclusively  concerned  with  "Molodensky's  problem,  "  the  gravimetric 
determination  of  the  earth's  physical  surface  and  external  gravity  field  from 
free-air  anomalies.  The  free-air  anomaly  Ag  is  defined  as  the  difference 
between  gravity  g  measured  at  ground  and  normal  gravity  y  referred  to  the 
telluroid.  The  discussion  of  the  solutions  of  Molodensky's  problem  is  in 
general  limited  to  "linear  solutions,  "  in  which  second  and  higher  powers  of  the 
elevation  (which  is  small  as  compared  to  the  dimensions  of  the  earth)  and  of 
the  terrain  inclination  are  systematically  neglected.  This  limitation  to  a 
"linear  approximation"  considerably  simplifies  the  comparative  study  of  the 
different  solutions  and  is  also  justified  for  practical  reasons. 

In  Molodensky's  problem  the  free-air  anomalies  are  assumed  to  be  given 
at  every  point  of  the  earth's  surface.  This  distinguishes  it  from  "Bjerhammar's 
problem,  "  the  determination  of  the  earth's  physical  surface  and  external  gravity 
field  from  discrete  gravity  measurements,  which  is  beyond  the  scope  of  the  present 


report.  We  have  also  disregarded  methods,  mainly  proposed  by  N.  K.  Migal, 


that  avoid  the  use  of  a  normal  gravity  field:  they  are  interesting  conceptually, 
but  less  convenient  practically  and,  as  shown  by  Monin  (1962),  essentially 
equivalent  to  the  usual  methods,  which  do  use  a  normal  field. 

Part  A  considers  various  integral  equations  for  the  present  problem, 
and  their  linear  solutions.  Much  of  its  subject  matter  can  be  found  in  the 
literature:  it  has  been  collected,  supplemented,  and  presented  from  a  unified 
point  of  view. 

Part  B  gives  a  deduction  of  the  linear  solutions  from  one  particularly 
simple  and  obvious  solution.  Here  I  have  attempted  to  show  that  all  linear 
solutions  that  I  was  able  to  find  in  the  literature,  and  some  more,  can  be 
derived  in  an  elementary  way,  in  the  sense  that  no  integral  equations  are 
needed  at  all.  In  this  way  the  relation  between  the  various,  apparently  so 
different,  solutions  is  clarified.  This  comprehensive  deduction  may  claim 
some  originality,  although  the  general  technique  was  anticipated,  for  a  particular 
case,  by  Arnold  and  indirectly  also  by  Molodensky. 

Part  C  considers  the  application  of  the  various  methods  of  solution  to 
gravity -dependent  quantities  that  are  needed  in  geodesy.  A  comparison  and 
evaluation,  both  theoretically  and  with  respect  to  practical  application,  is 
attempted. 

Mathematical  details  that  would  unduly  interrupt  the  presentation  have 
been  relegated  to  two  Appendices. 
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Formulas  that  might  be  new  are,  for  instance,  (63),  (106),  (261),  (274), 
(277),  and  (297b).  It  is  quite  possible,  however,  that  some  of  them  have  been 
found  before. 

The  three  parts  of  the  report  can  be  read  fairly  independently.  Readers 
interested  only  in  a  general  view  and  in  practical  application  may  limit  them¬ 
selves  to  Part  C;  and  those  who  are  afraid  of  integral  equations  may  start  with 
Part  B,  although  they  will  probably  have  to  refer  to  sections  1  and  3  for  a  full 
understanding. 
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A.  THE  INTEGRAL  EQUATIONS 


1.  SPHERICAL  FORMULAS 

For  later  reference  we  first  state  certain  well-known  integral  formulas, 
which  are  valid  if  the  boundary  surface  is  a  sphere. 

SPHERICAL  HARMONICS  EQUIVALENTS.  Let  r  (radius  vector),  0  (polar 
distance)  and  X  (longitude)  be  spherical  polar  coordinates  in  space.  Then  r  =  R  = 
6371  km  is  the  radius  of  the  mean  terrestrial  sphere,  considered  as  the  boundary 
surface.  Consider  on  this  sphere  a  surface  layer  of  density  X(0,  \)/2nk,  where  k 
is  the  gravitational  constant,  and  its  potential  Y(0,  X)  on  the  sphere.  These  two 
functions  are  connected  by  the  formula 


..  R  x  , 

Y  =  5  r-r  r  d0  • 


(i) 


where  Y  refers  to  a  fixed  point  P  and  X  to  a  variable  point  P'  that  carries  the 
surface  element  R2da,  a  being  the  concentric  unit  sphere;  l0  is  the  spatial 
distance  and  the  spherical  distance  between  P  and  P' ,  so  that 


L 3  =  2Rsin 


•  i> 

in 


(2) 


(Fig.  1). 

Representing  the  functions  X  and  Y  as  series  of  spherical  harmonics, 
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(3) 


X(8,  X)  =  £  X.  (9,  X), 

ISO 

Y(0,  X)  =  £  YB(ft  X), 

n=o 

equation  (1)  is  equivalent  to  the  relation 

2R 

Y»(0,  X)  =  XB(e,  X)  (4) 

between  the  harmonic  terms  of  degree  n. 

Similarly  Stokes' formula 

Y  -  s  If u  (5> 

a 

(Y  represents  the  anomalous  potential  and  U  the  gravity  anomaly)  is  equivalent  to 

Yn(0,  X)  =  ^  14(0,  X)  .  (6) 

To  obtain  the  formulas  inverse  to  (1)  and  (5),  we  consider  the  elementary 
inversion  of  the  corresponding  relations  (4)  and  (6): 

RXn  =(n+£)Yn  , 

RUn  =  (n  -  1)  Yn 


or 

(RX  -£Y4  =  nYB  , 

(RU  + Y4  =  n  Yn  . 

Both  equations  are  of  the  form 
Z,  =  n  Yn  , 

which  is  the  spherical  harmonics  equivalent  of  the  integral  formula 


(7) 


(8) 
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Z 


R^  rr  Y  -  YP 
2 ff  ^  £04 


da 


(9) 


(Heiskanen  and  Moritz,  1967,  sec.  1-18).  Outside  the  integral,  functions  (such 
as  Z)  always  refer  to  the  fixed  point  P  under  consideration;  inside  the  integral, 
functions  (such  as  Y)  are  in  general  variables  of  integration,  and  therefore  the 
value  of  this  function  at  P  is  then  explicitly  denoted  by  the  subscript  P  (such  as 


Yp). 


Now  the  inversion  of  (1)  and  (5)  is  straightforward.  By  means  of  (7),  (8), 
and  (9)  we  obtain 

x  -  jl  .  £  rr  Y ' Y’ 

A  —  —  JJ  ”7T' 


2R  2ir 


u.-jf-ir/ 


da  , 


(10) 


Y-  Yp 


da  . 


(ID 


APPLICATION  TO  THE  ANOMALOUS  GRAVITY  FIELD.  The  spherical- 
harmonic  expansion  of  the  anomalous  potential  in  space  may  be  written 

'O  ,  n  Nn  +  1 

T  =  £  (  T»  • 

v  r  - 


(12) 


first  and  second  degree  terms  being  omitted  as  usually.  Differentiating  we  find 


JB  +  1 


dT  00  RB 

-3T-.E  <n  +  1J  7**  T*  • 


The  gravity  anomaly  in  space  is  given  by 


dT  2T  ®  Rn  +  1 


We  write  this  as 


see,  e.g.,  (Moritz,  1964).  Its  solution  is  Stokes' formula  (5): 


T  -  If  AgS#)<to. 

a 

Another  integral  equation  is  obtained  by  representing  T  as  the  potential 
of  a  surface  layer.  Denote  the  product  of  surface  density  and  gravitational 
constant  by  cp;  then 


T  -  R2  f  J  f  da  . 


Inserting  this  into  the  boundary  condition 

dT  2T 
dr  R  ^ 


and  taking  account  of  the  discontinuity  of  the  normal  derivative  of  a  surface  potential 
we  obtain  tile  integral  equation  for  <p  , 

3R  nr <  <p  ,  1 

0  =  (21) 
see,  e.g.,  (Heiskanen  and  Moritz,  19o7,  sec.  8-b).  Inserting  (19)  into  (21)  we 


*-s  (*♦£)• 


Expressing  T  by  Stokes'  formula  (18)  we  obtain 


*>  *  Sr 


as  a  solution  of  the  integral  equation  (21). 
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GREEN'S  SURFACE  IDENTITIES.  Let  U  and  V  be  two  functions  defined 


on  an  arbitrary  surface  S.  Then  the  relation 


JJ  DM,  V)  dS  -  -  Ju 

Si  c 


av 

dl/ 


ds  -  JJ  UAe  VdS 
Si 


(24) 


holds.  Here  C  is  a  closed  curve  on  the  surface,  Si  is  the  part  of  the  surface  S 
that  is  enclosed  by  C,  ds  is  the  line  element  of  C,  dS  is  the  surface  element  of 
S,  v  is  the  tangent  to  S  that  is  normal  to  C,  D  is  Beltrami’s  mixed  differential 
parameter,  and  A2  is  Beltrami's  second  differential  parameter  for  the  surface. 

See  (McConnell,  1931),  p.  189,  eq.  (69);  the  notation  is  slightly  different. 

If  the  surface  is  referred  to  orthogonal  parameters  and  u2 ,  so  that 
the  line  element  is  represented  by 

ds2  =  hi2  dui2  +  h t  du22  ,  (25) 

then  the  differential  parameters  are 


D(U,  V) 

1 

-  -rs 

hi 

dU  dV  |  1 

dui  3ui  hs 

au  ay 

aus  aus  * 

(26) 

As  (F)  = 

1 

hi  he 

"  J-(  k  &F  >v 
_  dui  v  hi  3ui  J 

+  a  /'hi  aF-\ 
aueVhs  au,  /  ’ 

(27) 

Formulas  for  general  non-orthogonai  surface  parameters  are  given  in  (McConnell, 
1931)  on  p.  187. 

If  the  surface  were  a  plane,  then  ui  =x  and  u3  =  y,  and 


D(U,  V) 


dU  dV  dU  dV 

dx  dx  dy  by  ’ 


A2(F)  = 


(28) 

(29) 
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Hence  we  see  that  (24)  is  the  two-dimensional  analogue  of  Green's  first 
identity  in  space;  see,  e.g.,  (Moritz,  1964,  sec.  2.1).  Therefore  (24)  is  called 
Green's  theorem  for  a  surface. 

If  S  is  a  closed  surface,  and  if  the  curve  C  is  contracted  to  a  point,  then 
Si,  taken  to  be  the  exterior  of  C,  becomes  the  whole  surface  S,  and  the  first 
integral  on  the  right-hand  side  of  (24)  reduces  to  zero.  As  a  limit  there  remains 

fj  D(U,  V)  dS  =  -  JfUAaVdS  .  (30) 

S  S 

Because  of  the  symmetry  of  D,  also 

JJ  D(U,  V)dS  =  -  JJVAaUdS.  (31) 

S  S 


In  these  identities  the  functions  F,  U,  V  are  supposed  to  be  continuous  and 
twice  differentiable  in  the  region  considered:  Si  in  (24)  and  S  in  (30)  and  (31). 

The  function 

V=J,  (32) 

where  A  is  the  distance  of  a  variable  point  P'  from  a  fixed  point  P  on  the  surface, 
satisfies  this  requirement  everywhere  except  at  P,  where  l/l  becomes  infinite. 

In  this  case  it  may  be  shown  that  (31)  still  holds: 


whereas  (30)  must  be  slightly  modified: 

j;D(F,i)dS  =  -J'RF-F,)A3(i)dS.  (34) 

s  s 


II 

I 


As  we  shall  see  by  equation  (39),  A,  (1/1)  has  a  strong  singularity  (like  1/1=), 
Which  is,  so  to  speak,  neutralized  by  subtracting  from  the  function  U  its  value 
U,  at  the  fixed  point  P.  To  prove  (34)  and  (35),  exclude  first  P  by  a  small 
circle  C  of  radius  c ,  apply  (24)  and  let  subsequently  c  -  0. 

These  formulas  hold  for  an  arbitrary  closed  surface  that  is  sufficiently 
smooth.  Now  we  shall  specialize  these  general  formulas  to  the  sphere  r  =  R, 
whose  line  element  is  given  by 

ds2  =  R2  (dtf  +  sin2  0  d  A2 )  , 

so  that 

ui  =  0  .  u2  =  \  , 

R.  h2  =  RsinO  . 

Then  (26)  and  (27)  become 


(35) 


D(U,  V)  -  rAr  (  (I  V  f  —4~  U  V  ^ 

R  0  0  sin20  7a/’ 

As (F)  -  i  (  F  cot8  +  F„  +  *5  F  'l 

R  v  s  ee  siirQ  u  J  ■ 

where  U0  =  5U/5S,  U9f)  =  Da  U/56= ,  etc. 

Alternatively  we  may  refer  the  sphere  to  coordinates  *  (spherical  distance 
from  P)  and  or  (Azimuth  in  P);  see  Fig.  2.  These  coordinates  *  and  a  exactly 
correspond  to  8  and  A,  the  origin  being  now  the  fixed  point  P  instead  of  the 
north  pole.  Hence,  alternatively, 

D(U,  V)  =  -5  C  u  V  +— U—  u  V 

R  \  ip  $  sin‘0  a  a/’ 

^(F)  =  A  (  F  cor^  +  F.,  -f— 4  F  ^ 

R  V  ^  00  sin  0  roa  J  * 


(36) 
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For  the  sphere  we  have  t  =  Z0  =  2R  sin  ^ ,  so  that 


—  =  - V 

t0  2Rsin^  * 

(37) 

“'Or) --*£#£• 

(38) 

A2(*9  )  f83+4 FIT  * 

(39) 

This  follows  by  straightforward  evaluation  of  (36)  for  the  particular  function  (37). 
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2.  INTEGRAL  EQUATIONS  FOR  THE  ANOMALOUS  POTENTIAL 


In  (Moritz,  1964)  we  have  derived  two  different  but  equivalent  integral 
equations,  which  are  both  due  to  Molodensky.  We  shall  recall  the  n.ain  steps  of 
the  derivation,  which  we  also  need  to  derive  a  third  alternative  integral  equation. 
Our  starting  point  is  equation  (49)  of  (Moritz,  19b4): 


T±(i)-!!!-|dS  .  o. 

nil  \  l  S  i  }n  J 


(40) 


S  is  any  known  surface  approximating  the  physical  surface  of  the  earth,  such  as  the 
norn.al  surface  as  defined  in  (Moritz,  19o4,  sec.  2.3)  or  Hirvonen's  telluroid.  (In 
the  previous  report  cited  we  have  used  the  symbol  Z  instead  of  S,  but  this  notation 
is  son.ewhat  awkward.)  The  outward  directed  norn.al  to  S  is  denoted  by  n,  and  l 
is  the  spatial  distance  between  a  fixed  point  P  and  the  variable  surface  element  dS 
(Fig.  3). 

Consider  an  arbitrary  function  F  in  space.  Its  norn.al  derivative  may  be 
expressed  in  two  alternative  ways;  see  equations  (18')  and  (21)  of  (Moritz,  1964): 


dF  dF 

—  =  COS  P  --  -  D(F,  h)  cos  0  , 


9F  "  1  -  D(F,  h)  cos  p  , 


(41) 


(42) 


Sn  cos  P  dh 

where  h  is  the  elevation  of  the  terrain  and  is  its  angle  of  inclination. 

In  a  local  cartesian  coordinate  system,  the  z-axis  pointing  vertically  upward, 
the  x-axis  pointing  north  and  the  y-axis  pointing  west,  the  expression  D(F,  h)  is 
defined  as 
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*  ; 


D(F,  h)  =  —  —  +  2il 
3x  3x  ay  ay  ‘ 

The  expression  D(F,  h)  is  defined  in  the  same  way,  the  derivatives  along  the 
surface 


^f  _  aF  +  aF  ah  ^f  _  aF  aF  ah 

dx  ax  ah  ax  ay  ay  aii  ay 


(43) 


(44) 


replacing  the  horizontal  derivatives  c^x  and  ^y.  Other  expressions  for  D 
and  D  will  be  given  below.  The  vertical  derivative  of  F  is 


3F  _  aF 
3h  az 


(45) 


Equation  (41)  thus  expresses  the  derivative  of  F  along  the  surface  normal  n  (which 

is  not  in  general  vertical)  as  a  linear  combination  of  the  vertical  and  horizontal 
derivatives . 


For  F  =  T  we  have 


&  .  -Ag-4&T 

dh  Ag  y  5h  T’ 


ax 

ST 

dy 


y  ax 


t  =  . 


=  -yr\ 


(46) 

(47a) 

(47b) 


where  y  is  normal  gravity  and  £  and  r)  are  the  computers  of  the  deflection  of  the 
vertical.  Hence 


D(T,  h)  -  y  (£  tan  ft  +77  tan  0a  ) , 


(48) 


where 


3h 


tan^  tan/3a 


(49) 
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are  the  inclinations  of  a  north-south  profile  and  an  east-west  profile,  respectively. 
Inserting  (41),  with  D(T,  h)  expressed  by  (48),  into  (40)  we  obtain 


[±( 

'In  1  cosfi  I 

L  3n' 

^  I  y "  y  ah  i  J 

TdS  = 


2^  -  y  (£tan&  +T?tan02)]  cos£dS  , 


(50) 


which  is  our  basic  integral  equation  in  the  first  form;  it  is  identical  with  equation 
(48a)  of  (Moritz,  1964). 

Secondly  we  use  (42),  again  setting  F  =  T.  Molodensky  (Molodenskii  et 
al. ,  1962,  p.  85)  proved  that  for  two  functions  U,  V 


!T  D(U,  V)cos0dS  =  -  jr  U  A  a  V  cos  /3  dS 

s  s 

=  -  pr  VAa  Ucos  0dS  , 

j  u 

s 

where,  in  our  local  coordinate  system,  the  expression  A2F  is  approximately 
defined  by 

_  S32F  ^  Sa2F 


(51a) 

(51b) 


AaF  =  + 

3x2  df 


more  rigorous  expressions  will  be  given  below.  By  (42)  and  (46)  we  obtain 

=  — 1—( ^-Ag  +  -^T)-7  5(T,  h)cos  0  . 

I  3n  lcosB\  6  y3h  /  1  r 


ij3  v  “  y 

From  the  definition  of  D  it  follows  that 

-  'T 


}5(T,  h)  =  5({,h)-TD(j,») 


(52) 


(53) 


(54) 


Now  we  substitute  (53)  and  (54)  in  (40)  and  apply  (51a),  with  U  =  T /l.  The  result  is 
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(55) 


_  1  »(.  r  d  /'  1  \  1  Sy  1  —f  1  ,  \  . 

T  '  2n  g-  L  dn  ^  i)  '  y  bh  £cos/3  ‘  D(^  £  ’  h)cos^  ' 

-  A2  h  1  TdS  =  ^-  ff  — dS  . 

£  J  2 it  J  < J  £  cos  j3 


which  is  the  second  form  of  the  basic  integral  equation;  it  is  equation  (48b)  of 
(Moritz,  1964). 

A  third  form  is  obtained  by  using,  instead  of  (54),  a  different  transformation: 


7  D(T,  h)  =  -D(T,  h  -  hP) 

£  t 

*  6(t.  5-^)-(h-h,)D(i,  t) 

-  d(t,  5[(h  -h,)T,  j]+Td(|.  h).  (56) 

The  possibility  of  these  manipulations  follows  again  from  the  definition  of  D .  In 
agreement  with  (51)  we  have 

JjD(T,~-^)cos/3dS  =  -  A2Tcosj3dS  ,  (57) 

S  S 

JJ  D  j~(h  -  hp ) T,  7  I  cos /3dS  =  -  JJ  (h  -  hP)TA2  (  7  )cos  j3dS  . 

S  S  (58) 


Again  we  substitute  (53)  in  (40),  but  now  we  apply  the  transformation  (56)  and  use 
(57)  and  (58).  The  result  is 

y  dh  Tzkj  +  (h~h')A3{^i)cosP 

J 

+  D(i.  h)cos^  J  TdS  =  £  yAgjJh^A^Wjds. 


+ 


£  cos/3 


(59) 
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This  is  a  third  form  of  the  basic  integral  equation. 


Alternative  expressions  may  be  obtained  by  applying  (42)  to  the  function 


9  r  l  >  l  hr  l  n  i  n 

— (  7  i  - -  rr  (  7  )  -  D{  -  ,  h  )  cos  0  . 

9n^£y  cos/3  3h  V  i  /  vr  s 

On  this  substitution  the  basic  equations  (50),  (55),  and  (59)  become 

T  -  ~  rr  r  — ^  ~  (  7  )  -  -  '  d(  7  ,  h ")  cos/3  1  TdS  = 

2tt  -h  l  cos  8  9h  V  £  y  y  9h  £  v  £  J  J 


=  ^  rr  |  [  Ag  -  y  (^tan/3x  +r)tan/3s)]  cos/3dS  , 


.  JL  rpf _ L_  ±(k\„i*2  — L 

2rr  L  cos/3  9hV£y  y  9h  £cos 


-  -25(},h)cosS 


^  i2h]  tus  =  i  rr  as. 


2rr  £  cos  /3 


T  -  J-  r  -L-  rj  V  -  ■■  +(h  -  hp)Ar  ^  7  )  cos/3]t  dS  = 

2tt  -  *  l  cos  B  9h  s  £/  y  9h  £  cos  /3  .  £  J 


»  ’.'T^f^B-Orl'OA^cos^ldS. 

s 


These  equations  are  more  explicit  and  therefore  better  suited  for  practical 
solution;  the  third  form  has  even  become  simpler  because  the  D  term  cancels  out 
in  (63). 

The  disturbing  potential  T  may  be  determined  from  any  of  these  integral 
equations.  The  necessary  data  are:  Ag  and  the  deflection  components  £  and  77 
in  (61);  Ag  only  in  (62);  and  Ag  and  A2T  (which  is  essentially  the  anomalous 
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gradient  of  gravity)  in  (63).  Equations  (61)  and  (62)  are  due  to  Molodensky;  (63) 
might  be  new. 

The  solution  of  these  equations,  including  practical  aspects,  will  be  discussed 
in  sec.  4. 

THE  OPERATORS  D,  Dt  AND  A2  .  First  of  all,  the  reader  should  be 
warned  that  the  operators  D  and  A2  as  used  in  this  section  are  not  absolutely 
identical  with  those  denoted  by  the  same  symbols  in  sec.  1,  although  they  are 
closely  related.  Therefore,  the  integral  formulas  (30)  and  (51)  are  very  similar 
but  not  identical;  note  the  factor  cos  fi  in  (51). 

The  reason  is  that  in  the  present  section  the  differential  parameters  D,  D, 
and  A 2  (even  when  referring  to  a  surface  S)  are  essentially  defined  in  terms  of  a 
three-dim en sional  coordinate  system.  It  is  convenient  to  use  some  simple  system 
of  orthogonal  curvilinear  coordinates  qx,  q 2,  q3  ,  such  that  the  surfaces  q3  =  const, 
are  approximately  level  surfaces,  qx  is  approximately  the  latitude,  q2  is  approxi¬ 
mately  the  longitude  and  q^  is  approximately  the  elevation.  Tlien  the  line  element 
will  be  given  by 

ds2  =  hi 2  ckji 2  +  l*,2  dq22  +  h32  dq32  ,  (64) 

where  hi ,  h2,  1^  are  functions  of  qx ,  q2 ,  qa  (they  should  not  be  confused  with  the 
elevation  h).  Then  the  following  definitions  hold  (Molodenskii  et  al. ,  1962, 


pp.  83-85): 


D(u-  v>  - 


5V  l  ha  5U  dV 
hg2  Sq2  dqs  ’ 


(65) 


D(U,  V) 


=  Ji3S2USaV+JbSaU&V 

hr  aqi  fch  9Q2  aq2  ’ 


(66) 
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(07) 


a  (F)  =  JL_ 

2 1  '  ^  ^  L  Sqi  V  h:  dqx  J  V  hg  dog  J  J  ' 

As  we  have  seen,  dg  F  denotes  differentiation  along  the  surface.  If  F  is  a  function 
defined  on  the  surface  S  only,  F  =  F(qi ,  qg),  then 


M  _  ^F(qx ,  (fa)  . 

aqi 


(68a) 


if,  however,  F  is  originally  a  function  defined  in  space,  F  =  F(qx ,  qa ,  qa),  and 
if  the  surface  S  has  the  equation  qa  =  h^  ,  qa ),  then 
2a£  =  bF(ch  ,  <b  ,  hfa  ,  qa  )  ) 

3F(0l ,  Qa  ,  q3)  +  ^F(q:,  Qs ,  qa)  dh 


or  briefly, 


aaF  _  ^F  t  dF 
dqj  dqj  3q:.  5qi 


5q3 


(68b) 


For  surface  functions  such  as  the  elevation  h  or  the  reciprocal  distance  1  /£,  defini¬ 
tion  (68a)  holds;  for  space  functions  such  as  the  anomalous  potential  T,  definition 
(68b)  n  ust  no  applied. 

The  definitions  (66)  and  (67)  reduce  to  (26)  and  (27)  only  when  the  surface  S 
under  consideration  is  a  surface  qa  =  const,  (with  ha  =  l,  which  can  always  be 
achieved  bv  a  suitable  choice  of  coordinates  q-s ).  Then  $  =  0  and  (51)  reduces  to 
(30)  or  (31).  Broadly  speaking,  this  holds  when  S  is  approximately  a  level  surface, 
which  the  physical  surface  of  the  earth  obviously  is  not. 

The  conditions  to  oe  satisfied  by  the  functions  entering  in  (51)  are  the  same 
as  for  (30)  and  (31):  the  functions  U  and  V  are  supposed  to  be  continuous  and 
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twice  differentiable  over  the  whole  surface  S.  The  function  F  =  1/f  has  a 
singularity  at  P  and  does  not  therefore  everywhere  satisfy  these  conditions. 

In  spite  of  this  fact,  eq.  (51a)  still  holds  for  U  =  1/f,  whereas  (51b)  must  be 
modified  by  replacing  V  by  V  -  V„ ,  in  order  to  neutralize  the  strong  singularity 
of  As  (1/f).  This  is  in  complete  analogy  to  (33)  and  (34)  and  is  proved  in  the 
same  way.  This  fact  was  used  in  the  transformations  necessary  to  derive 
(55)  and  (59). 

As  an  example,  we  consider  the  geodetic  coordinate  system  <p  (ellipsoidal 
geographic  latitude),  X  (ellipsoidal  longitude),  and  h  (height  above  the  reference 
ellipsoid).  The  line  element  in  these  coordinates  is  given  by 


ds2  =  (M +  h)2  c*p?  +  (N +  h)2  cos2  odA2  +dh2 


(Molodenskii  et  al. , 

M  = 


1962,  p.  10),  where 
a2  b2 

(a2  cos2  (p  +  b2  sins  o  )3^ 


(69) 


(70) 


is  the  meridional  radius  of  curvature  and 


N  = 


(a2  cos2  +b2  sin2 


(71) 


is  the  east-west  radius  of  curvature  of  the  ellipsoid,  whose  semi-axes  are  a 
and  b.  Hence 


Qi  =  <P  ,  Q3~h; 

hx  =  M  +  h,  h2  =  (N+h)cos<p,  h3  =  1  . 
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Equation  (66)  thus  becomes 


on,  V)  =  — L_,  ^  iay  M 

'  *  (M+hf  dp  ap  (N  +  h)P  cos3 p  dX  dX 


(72) 


For  (67)  we  find 


(73) 


Aa(F)  =  - i - ris.r (N+h)cos<p  3*F\  + 

(M  +  h)  (N  +  h)cosp  L dp V  M+h  dp  J 

+  ia  r _ M  +  h  . 

dX  v.  (N  +  h)cos<p  dX  J  \  ’ 

Differentiation  of  M  and  N  with  respect  to  p  is  straightforward  (their  derivatives 
with  respect  to  X  are  even  zero),  but  it  should  be  noted  that  according  to  (68b) 


ik  ( 

'(N  +  h)cosp  N 

_  _d_ 

f  (N  +  h)cosp' 

ui  ( 

dp  ' 

^  M  +  h  J 

dp 

V  M+h  „ 

)  dhV 

Jk  ( 

'  M  +  h  \ 

d 

/■  M+h 

\  dh 

dX^ 

^  (N+h)  cosp 

dh 

k  (N  +h)  cos  p  > 

'  dX  * 

(74) 


If  the  reference  ellipsoid  is  considered  as  a  sphere  (spherical  approximation 
see  next  section),  then 

M  =  N  =  R  , 
and  we  are  left  with 

1  f'  dg  u  dsjV  A  1  dfeU  d2  V 


D(U,  V)  =  ^  ^  21  +  -4-  ^ 

K  v  -.p  dp  cos  p  dX  dX  y 


As  (F)  =  j?  (j  —■  ten®  +  2^  + - 4 M  )  , 

K  V  dn  d«  cos  m  d  A  y 


(75) 

(76) 


dp  dp"5  cosdp  dXJ 

in  formal  analogy  to  (35);  note  that  for  the  sphere 

p  =90 P  -  0  .  (77) 

In  the  local  cartesian  coordinate  system  we  have  the  expressions  (43),  (44) 
and  (52)  already  given. 
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3.  VARIOUS  APPROXIMATIONS 


SPHERICAL  APPROXIMATION.  The  geodetic  integral  equations  become 
simpler  and  can  be  solved  more  easily  by  applying  certain  approximations  which 
are  permissible  from  ihe  point  of  view  of  accuracy. 

Usually  i  the  gravity  field  of  an  ellipsoid  of  revolution  is  taken  as  the 
normal  field.  The  flattening  of  any  suitable  reference  field  is  small.  Hence, 
when  dealing  widi  ellipsoidal  quantities,  it  is  convenient  to  use  series  expansions 
with  respect  to  the  flattening  f  or  a  similar  small  parameter. 

Our  integral  equations  deal  with  quantities  of  the  order  of  T,  which  are  very 
small  themselves.  In  these  integral  equations  it  is  therefore  possible  to  neglect 
all  terms  containing  f,  f ,  etc.,  as  a  factor  of  T.  This  is  the  spherical  approxi¬ 
mation.  The  error  in  the  height  anomalies  £  (sec.  15)  of  this  spherical  approxi¬ 
mation  is  thus  of  the  order  of 
f £  »  0.003£  . 

If  £  =  100  m,  this  amounts  to  0. 3  m. 

A  convenient  visualization  of  the  spherical  approximation  is  furnished  by 
plotting  the  heights  h  above  the  reference  ellipsoid  as  the  heights  above  a  mean 
sphere  of  radius  R;  see  Fig.  4.  It  should  be  kept  in  mind,  however,  that  this  is 
only  a  visualization,  which  has  no  counterpart  in  reality .  It  would  be  plainly  wrong 
to  say  that  the  earth  is  referred  to  a  sphere  instead  of  an  ellipsoid,  at  least  in  the 
literal  geometrical  meaning. 
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Consider  a  definite  case.  The  true  radius  vector  of  P'  (Fig.  4a)  is  given 


by 


rt  r  u  e  -  a  (1  -  i  e  3  sin  <p  +  •  •  •  )  +  h 

=  R(l  +  Te/2'£e/2  sin2<p  +  •••)  +  h 


(78) 


where 


_  /a  _  a2  -  b2 

is  the  square  of  the  second  excentricity  and 

R  =/?b 

is  the  mean  radius  of  the  earth.  Neglecting  e/2  we  get 
r  =  R  +  h 

as  the  spherical  approximation  of  (78).  This  corresponds  to  the  geometrical 
interpretation  of  Fig.  4b. 

In  agreement  with  (78  ),  we  have  to  the  same  approximation 

i_= 

dh  or 


(78') 


(79) 


and 


with 


t  -  (r2  +rp3  -  2rp  rcos 0  )~^ 


r  =  R  +  h  , 
rP  =  R  +  hp  . 


(80) 


(81a) 

(81b) 


27 


Hence  we  btain  by  differentiating  (80) 

r -  r ■  rp cos0 
ahvly  '  br\is  '  "  p  * 

5  /"  _1  \  _  5  /"  _1  "N  _  rp  -  rcos$ 

aiv  \ i  J  arP  \  i/  jj3 


This  is  easily  transformed  into 

1/ JA  =  .  J_  +rp8  -  f 

ahv£y  2r  l  2rP  * 

1^  — 


ah' 

a 

bhp 


Cl) 


2rP>:  2rP  l 3 

From  the  spherical  approximation 


(82a) 

(82b) 


kM 

y  =  -pr 

we  further  find 

i  il  =  .2 

y  ah  r  * 

Hence,  in  (62)  and  (63)  we  may  set 


(83) 


a  /  1  ^  1  5y  1  _  3  |  rps  -  r? 

bh  v  £ )  y  dh  l  2r£  2r£3 

The  spherical  approximation  of  the  operators  D  and  As 


(84) 

is,  of  course,  given  by 


(75)  and  (76). 

PLANAR  APPROXIMATION.  Comparing  (78)  and  (78 ')  one  might  object  that 


it  is  illogical  to  neglect  the  term 

-  5  sin*  e,z  R  , 
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which  may  become  as  large  as 

-  \  e'2  R  =  -  14  km  , 

o 

when  at  the  same  time  h,  which  is  8  km  at  its  maximum,  is  retained. 

That  this  objection  is  not  valid  may  be  seen  by  considering  the  consequences 
of  completely  neglecting  h  .  Then  the  earth's  surface  would  be  represented  as  a 
sphere,  and  all  mountains  and  hills  would  be  leveled.  This  would  amount  to 
neglecting  also  the  inclination  of  the  terrain,  which  may  attain  45'  and  more 
in  steep  hills  and  mountains.  This  is  clearly  inadmissible. 

So,  while  it  would  be  permissible  to  neglect  the  elevation  as  such,  it  is  not 
possible  to  do  so  on  account  of  the  rapid  changes  of  h,  which  cause  the  inclination. 
In  other  words,  one  may  neglect  h  whenever  it  does  not  enter  through  its  horizontal 
derivatives 


dh 


=  tan  ft  , 


dh 


=  tan/32 


dx  r'1  '  dy 

or  through  similar  expressions. 

Such  an  expression  is 
h  -  h„ 

K  ’ 

where  l0  is  the  chord,  corresponding  to  l  ,  of  the  sphere  r  =  R  (Fig.  4b): 
le  =  2Rsin|. 


(85) 


(86) 


For  small  distances  ,  the  quantity  (h  -  hp  )/£,  is  of  the  order  of  the  inclination 
/3,  as  Fig.  5  shows. 
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1 1) 

where  l0  =  2Rsin-^-.  This  is  the  spherical  approximation  of  (88).  By 
neglecting  h  except  in  the  term 


we  obtain 

l9  =  **|><^)  ]  =  io8+(h-hpf  (88') 

as  the  planar  approximation  of  (88). 

As  a  second  case,  consider  (84),  which  may  also  be  written 

_d_  /^.l^  1  dy  1  3  (h  -  hf, )  (r  +  rp ) 

bh  \lJ~y  oh  £  3  2r£  2r?  * 

The  planar  approximation  is  obtained  by  putting  r  =  rP  =  R: 

1 V  1  ^  —  h_lV 

ahV^y’yahc  2R i~  p  »  (89) 

where  l.  is  expressed  in  terms  of  l0  by  (88'). 

Neglecting  h/R  within  the  parentheses  of  (87)  may  be  interpreted,  in  a 
formal  mathematical  manner,  as  letting  R  -  ®  (within  the  parentheses!),  that  is, 
as  performing  the  formal  transition  from  a  sphere  to  a  plane.  The  expression  (85) 
will  Become  large  only  if  l0  is  small,  that  is,  in  the  neighborhood  of  P.  In  this 
neighborhood,  however,  the  sphere  may  be  replaced  by  its  tangent  plane.  Similarly, 
(88")  resembles  the  plane  formula  of  Pythagoras.  These  facts  indicate  that  the 
planar  approximation  may  be  visualized  in  the  way  of  Fig.  6:  the  elevations  h 
above  sea  level  or  above  the  ellipsoid  are  superimposed  on  a  plane.  This  inter¬ 
pretation  furnishes  a  convenient  name  for  the  approximation  considered. 
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This  name,  planar  approximation,  however,  should  not  be  taken  too  literally. 
The  transition  to  the  plane,  R  -®,  can  only  be  performed  in  expressions  of  the  order 
of  the  elevation  h,  but  nc  generally. 

We  may  thus  summarize  the  different  approximations  that  are  consistent 
from  the  point  of  view  of  accuracy: 

Normal  field  (normal  potential  U,  normal  gravity  y,  etc.):  strictly 
ellipsoidal  formulas. 

Anomalous  field  (anomalous  potential  T,  gravity  anomalies  Ag,  etc.): 
main  part  (not  depending  on  h):  spherical  approximation, 
correction  terms  (order  of  h):  planar  approximation. 

In  view  of  the  conceptual  difference  between  spherical  and  planar  approxi¬ 
mation  it  is  remarkable  that  the  accuracy  requirements  corresponding  to  the 
spherical  approximation  for  T,  etc. ,  entail  the  planar  approximation  for  the 
correction  terms  in  a  quite  conclusive  manner. 

LINEAR  APPROXIMATION.  Using  the  planar  approximation,  we  may 
expand  the  expressions  involved  into  power  series  with  respect  to 

—J  -  ,  tan/3  ,  (90) 

and  similar  quantities.  Neglecting  terms  of  second  and  higher  degree  in  these 
quantities  and  retaining  only  the  linear  terms  is  called  linear  approximation. 

The  linear  approximation  to  cos/3  is 

cos/3  =  1  .  (91) 
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because 

cos  0  =  (1  +  tan8  0)  ^  =  l-g  tan8  0  +'  *  * 
differs  from  1  only  by  terms  of  second  and  higher  degree  in  tan  0. 

As  a  linear  approximation,  we  further  have 

A  =  to  (92) 

for  the  same  reason,  according  to  (88"). 

Thirdly  we  use  the  surface  element  dS,  which  we  shall  also  need  subsequently, 
to  further  illustrate  the  meaning  of  the  various  approximations.  Fig.  7,  in  which 
the  profile  is  taken  along  the  maximum  inclination  0,  shows  that,  as  a  spherical 
approximation, 

dS  cos  0  =  r8  do  , 

where  do  is  the  element  of  solid  angle  (surface  element  of  unit  sphere).  Hence 
we  have: 


spherical  approximation: 

dS  =  r8  sec  0  da  , 

(93a) 

planar  approximation: 

dS  =  R8  sec  0  da  , 

(93b) 

linear  approximation: 

dS  =  R8  da  . 

(93c) 

The  linear  approximation  takes  account  of  the  main  part  of  the  effect  of 
topography.  Usually  it  should  be  sufficient  in  practice.  Furthermore,  the 
relation  between  the  different  solutions  as  considered  in  the  present  report  is 
most  obvious  when  the  discussion  is  limited  to  the  linear  approximation,  the 
relations  for  higher  approximations  becoming  rapidly  more  complicated.  There¬ 
fore,  we  shall  here  limit  ourselves  to  linear  solutions. 
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4.  LINEARIZATION  AND  SOLUTION  OF  THE  INTEGRAL  EQUATIONS  FOR  T 

By  means  of  the  spherical,  planar,  and  linear  approximation,  using 

equations  (89),  (91),  (92),  and  (93),  we  transform  the  basic  integral  equations 

(61),  (62),  and  (63)  into 

_  Ra  f.f'T  3  h  -  h*  -r  i 

T  "  2rr  ^  L  2R£0  '  /0  ’  h)]TdcJ  = 

Ra  i  i  r  i 

=  2^  jj* 7"laS  ~y  +r? tan£a)  |  do  , 

o  0 

Ra  -f  3  h  -  hp  r  1 

1  -  7T  1  1  I  * 


(94a) 


2rr  -  L  2R£n 


2D 


<  1  \  A2  h  i 

ut-v-t-]™- 


_  Ra  r>  r  Ag 

-2TT  ^  £a 


do  , 


(94b) 


Ra 


T-l?  J -r [ 2^r  •!li^+(h'h')AsCi')]Tda  = 


=  2p  If  7"  [Ag  -  (h  -  h,)  As  T]  da  . 


(94c) 


The  operators  D  and  Aa  in  these  equations  may  be  identified  with  those 
of  sec.  1,  defined  by  (26)  and  (27)  or,  as  a  spherical  approximation,  by  (35).  The 
reason  is  that  (35)  differs  from  (75)  and  (76)  only  by  d2F/d<p  being  replaced  by 
3F/o</),  etc..  Now,  for  surface  functions  such  as  l/l0  or  h,  these  two  derivatives 
are  identical,  according  to  (68a).  For  space  functions  such  as  T,  they  differ  by  a 
term  linear  in  h,  according  to  (68b),  which,  when  multiplied  by  h  -  hp ,  becomes 
quadratic  and  therefore  negligible  in  the  linear  approximations.  This  happens  for 
the  expression  (h  -  hp  )A3  T;  all  other  operators  D  and  Aa  involve  only  surface 
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•!  i 


1 


functions.  TTtfs  proves  our  assertion. 


Hence  we  may  use  (38)  and  (39).  The  integrand  on  the  right-hand  side  of 
(94c)  then  becomes 


3T 

2R£0  ’ 


because,  as  a  planar  approximation,  the  term  (h  -  hp  )/R  may  be  neglected. 


All  these  facts  being  taken  into  account,  the  integral  equations  (94a,  b,  c) 
become  finally 


t 


where 


Pi  =  •  y  (£tan0i  +  rjtan& ) , 

R2  « •»  1  f.  ,  .  dh  N  _  , 

“■•'a?  J-77^h'h>  's“*s)Tdo; 


Ps  =  ‘  T A2  h  , 


(97) 


R  *  **  1  /!  ,  oh  "n 

J:  -77-  (h  -  ",  -  2311^  —  Tdo  ; 

a 

p3  -  -  (h  -  hP  )A2  T  , 

Q3  =  0 

are  small  correction  terms  of  linear  order. 

They  are  readily  solved  in  two  steps.  As  a  first  step,  we  disregard  q, 
and  determine  an  approximation  T0  of  T  from  the  equation 

To  3R  rf^  da  =  ^  JJ  da  . 

4  IT  ^  *  1L  2ir  *  *  £ 


a  a 

This  is  a  purely  spherical  integral  equation  of  type  (17),  whose  solution  according 
to  (18)  is 


(98) 


T0  =~  JJ  (Ag  +Pl )  S  fy)da  . 


4  IT 

a 

We  then  put 

T  =  T0  +TX 

and  substitute  in  (96),  obtaining 


(99) 


(100) 


T-  -f  *  +  T,  -f  |  JJ  . 


o 


2tr  JJ  4 

a 
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TTie  first  two  terms  on  the  left-hand  side  cancel  with  the  first  term  on  the  right 
hand  side  because  of  (98),  and  there  remains 


3R  r»  p  T\ 

T'  II  71  *  A. 


4ir  JJ  4 
a 


(101) 


This  integral  equation  for  Tj,  is  of  type  (21),  so  that  its  solution  according  to  (23), 
replacing  Ag/2tr  by  qt,  is 


Tl  =  q,  JJ  q,  S  U»)  da  .  (102) 

o 

Combining  T9  and  Tx  we  have 

T  =  *rr  +  +  (103> 

o 

Consider  now  the  third  term  between  the  parenthesesof  (103),  3qt /2R. 

Since  q,  is  divided  by  R  in  this  term,  it  is  tempting  to  neglect  it  as  a  planar 
approximation  by  formally  letting  R  This  formal  argument  is  not  entirely 
convincing,  but  the  fact  that  3qt  /2R  is  negligible  nevertheless  holds  true;  the 
proof  is  rather  lenghty  and  will  be  deferred  to  Appendix  1 .  Hence  we  are  left 
with  the  simple  formula 

T  =  ~  j7(Ag  +  p,)S(0)do+q1  .  (104) 

O 

Using  (97)  we  may  as  follows  summarize  the  equivalent  linear  solutions  found: 


T  -  rj  [Ag  -  y  ( £  tan  &  +T7tan/32)]  S(0)do 

o 


(105) 
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(106) 


T  =  ^  JT  (Ag  -  TAS  h)  S  (0 )  da 
a 

'I  j4(li-h,  -  2sin*||)Td(,; 
a 

T  =  ^-  JJ  [Ag  -  (h  -  hp)AaT]  S(0)da  .  (107) 

a 

A  "zero -order"  approximation  to  T  is  given  by  Stokes'  integral 

jj?  j7  Ag  S  (^  )  do  ; 
o 

obviously  this  value  may  be  used  to  represent  i'  in  the  small  correction  terms 
Pj  and  qt . 

Expressions  for  second-order  terms  were  given  by  Arnold  (1959a)  and  Koch 

(1965). 
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5.  THE  INTEGRAL  EQUATION  FOR  THE  SURFACE  LAYER 
A  considerably  simpler  integral  equation  was  obtained  by  Molodensky  in 
an  indirect  way,  which  is  familiar  from  an  application  to  the  simpler  boundary- 
value  problems  of  potential  theory  (Kellogg,  1929). 

The  anomalous  potential  T  is  expressed  as  the  potential  of  a  surface  layer, 
or  coating,  on  the  earth's  physical  surface  S  (or  on  any  known  surface  close  to  it 
such  as  the  telluroid): 

T  =  JJ  \  dS  ,  (108) 

S 

where  <p  is  the  density  of  the  layer  (it  incorporates  the  gravitational  constant). 

It  is  known  from  the  theory  of  surface  layers  that  the  potential  and  its 
tangential  derivatives  are  continuous  on  the  surface  S,  whereas  the  derivative 
along  the  normal,  dT/dn,  is  discontinuous  there:  if  we  approach  S  from  the 
outside,  the  limit  of  oT/Sn  on  S  is 
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of  the  discontinuity  -  2ir<p  of  the  normal  derivative  and  the  continuity  of  the 
tangential  derivatives.  Let  this  direction  be  the  vertical;  then  /3  is  the  angle  of 
inclination  of  the  terrain.  Hence  the  limit  of  the  vertical  derivative  on  S  is 

(ID;  ’ 2  cos  0+  SET  ( l) dS '  (I10) 

u 

By  substituting  (108)  and  (110)  into  the  boundary’  condition  (46)  we  find 

2H.PCOS0- -(£  |f)  jJuxS  =  Ag.  (Ill) 

O  P 

which  is  to  be  considered  as  an  integral  equation  for  the  unknown  surface  density 
<p,  the  gravity  anomaly  Ag  at  the  earth's  surface  being  given. 

In  agreement  with  our  notational  convention,  quantities  outside  the  integral 
always  refer  to  the  fixed  point  P;  those  inside  the  integral  refer  to  the  variable 
surface  element  dS  unless  they  are  marked  by  the  subscript  P. 

If  we  compare  the  above  integral  equation  with  those  obtained  previously, 
that  is  (61),  (62),  and  (63),  we  see  that  (111)  is  much  simpler,  but  we  also 
recognize  an  essential  difference:  in  (111)  the  differentiation  d/5h  is  referred 
to  the  fixed  point  P,  whereas  in  the  former  equations  it  was  referred  to  the 
variable  point  dS. 

SPHERICAL  APPROXIMATION.  In  agreement  with  (84)  we  have 

S  ary  i  =  3  /-i*. 

dhfA-ey  \ybhsfi  2rP£  2rP£3  * 

because  l  is  symmetric  with  respect  to  r  and  rp .  Using  this  expression  and  (93a) 
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we  find 


y-  Q  J3  „  „  2 

2trpcos/3  -  - — 7  + - ^  ^  i*  sec/3<pdcr  -  Ag 

ua  ^2rPl  2rp;3  >/ 

as  the  spherical  approximation  to  (111). 

PLANAR  APPROXIMATION.  In  agreement  with  (89)  we  have 

9  ^  1^  r  1  I  =  _JL  +  h  ~  hp 

3hp  Vf/'  Vy  dh  ,£  2RX  £3  * 

and  the  surface  element  is  given  by  (93b).  Hence 

„  a  ,,r3R  ^  R2(h  -  hp) 

2ir<pcos0  -  JJ  1 _Yi  +  - 7a -  secptpdo  =  Ag 

cr  1  J 

is  the  planar  approximation  of  our  integral  equation. 

LINEAR  APPROX IMATION.  By  setting  i  =  £0 ,  cos/3  =  sec/3  = 
the  integral  equation 

3R  'v>  <P  ,  .  .  r»a  r> n  h  *  hp 

2inp  -  —  JJ  "p  cfcr  =  Ag  +  Rs  i  j  -j-T-  (pda 

a  0  o 


or 


where 


3R  |-»f>  (p  .  1  ,  n  v 

•  ‘  «  77  ^  -  ¥  <lig+Gl)  > 

cr 


Gi  =  r2  J  J  — 7~T  < 0  dcr  , 
ct 


and  (108)  becomes 

T  =  rS  j3‘  f  . 


(112) 


(113) 


1  we  obtain 


(114) 


(115) 


(116) 


(117) 
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where, as  usually,  L0  =  2Rsin(^/2). 

SOLUTION.  Comparing  (117)  with  (19)  and  (115)  with  (21)  we  see  that  the 
corresponding  equations  are  identical  if  Ag  is  substituted  by  Ag  +  Gx .  Making 
this  substitution  in  (18)  we  therefore  obtain 


T  =  —  r;  (Ag  +  G1)S(0)d<j 


(118) 


as  the  solution  of  the  system  of  equations  (115)  and  (117). 


Since  Gx  is  a  small  correction  term,  we  may  in  (116)  approximate  by  (22), 


so  that 


„  _  Ra  -  h  -  hP  f  3T  \ 

1  "  2 It  -  iaa  VAg  +2R  )d°  ' 

a 


(119) 


In  Appendix  I  it  is  shown  that  neglecting 


R5  r>  r>  h  ~  hp 


cr 


Cr) 


is  consistent  with  the  planar  approximation.  Hence  we  are  simply  left  with 


„  Ra  r'f'  h  -  hB  .  , 

Gl  =  J  J  Ag  **  ’ 


(120) 


Thus  the  combination  of  (118)  and  (120)  constitutes  another  linear  solution 
of  Molodensky's  problem. 

A  very  elegant  method  for  finding  higher-order  approximations,  which  is 
also  applicable  to  other  geodetic  integral  equations,  is  described  in  (Molodenskii  et 
al.,  1962,  pp.  120-3). 


6.  BROVAR’S  GENERALIZATION 


Molodensky  expressed  the  anomalous  potential  T  as  the  potential  of  a 
surface  layer  (108): 

T  =  ff  $  dS.  (12D 

S  1 

The  formal  reason  why  this  is  possible  is  that  1  /l  is  harmonic  as  a  function  of  P, 
and  therefore,  according  to  the  theory  of  linear  partial  differential  equations,  T 
is  also  harmonic  (that  is,  a  solution  of  Laplace's  partial  differential  equation). 

Brovar's  (1963c)  idea  is  to  replace  l/X  by  a  different  harmonic  function  E, 
arriving  at 


T  =  ff  <pEdS. 

s' 


(122) 


This  representation  is  valid  since  T  will  be  harmonic  if  E  is,  for  the  same 
reason  as  above.  We  may  consider  (122)  as  the  potential  of  a  generalized  surface 
layer,  and  the  surface  function  <p  as  a  generalized  surface  density. 

Then,  outside  the  surface  S  we  have  as  a  spherical  approximation 


5T  2T  -  r  3E 
g°  3rP  rP  dr„ 


25)  ds. 

rp  J 


(123) 


The  function  E  may  be  selected  in  such  a  way  that  the  "kernel” 


Y  =  .  dE  2E 
r0 

has  a  suitable  form. 


(124) 


For 


E 


4b 


(125) 


we  had,  according  to  (112), 


r 


r  -  r. 


1 3  > 


2r„  C  2rp  r 

that  is,  a  linear  combination  of  the  simpler  functions 


(126) 


1 


r2  -  r,2 


7  and  e 


(127) 


We  shall  now  try  to  find  functions  that  are  proportional  to  only  one  of  these 
functions  (127). 

To  arrive  at  a  suitable  generalization,  we  may  start  from  the  spherical- 
harmonic  expansion 


oo  rr' 

7  "  £  rrrrr  pr.  (cos  <J> )  . 

t  n  =  0  r? 


(128) 


considered  as  a  harmonic  function  of  the  point  P  in  space.  The  series  (128) 
remains  harmonic  when  the  individual  terms  are  multiplied  by  constant 
coefficients  \  ,  provided  it  converges.  Hence  we  may  use  the  representation 


E  =  S  \  P„  (eos  4>)  . 

n  =  0  rp 


(129) 


Then  the  kernel  (124)  becomes 


K  =  t  (n  -  1)  \  —its  PB  (cos  0 ) 


n=  0 


(130) 


According  to  (82b)  »ve  have 

-  t — = .  2  J-ri 


r,r 


2  J-(i').  JL 

3rp  V  SL  J  rP£ 


(131) 
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By  means  of  (128)  we  obtain 


r  ,r 


Hence  if  we  wish  to  have 


3  3 

r  "  7§~  =  E  (2n  +  1)  ~nr  Pu  (cos  $  )  . 


»=0 


(132) 


K  =  - 


a  z 
r  -  r. 


(133a) 


rp  l 


we  must  take 


= 


2n  -l-  1 
n  -  1 


(n  ^  1) ; 


for 


K  = 


rP  f 


(134a) 


(133b) 


we  have 


\  = 


1 

n  -  1 


(n  +  1) 


(134b) 


This  is  seen  by  comparing  (130)  with  (132) and  (128),  respectively. 

The  coefficient  \  remains  undetermined:  die  value  n  =  1  is  an  "eigenvalue. 
This  gives  rise  to  slight  complications,  which  will  be  skipped  over  here;  see 
(Brovar,  1964b)  and  also  (Moritz,  1965a). 

Both  cases  (134a,  b)  will  now  be  considered.  The  first  case  is  of  great 
practical  significance,  whereas  the  second  presents  several  features  of 
theoretical  interest. 

THE  KERNEL  (rc  -  r.s)/*a.  In  agreement  with  (134a)  we  take 
1  £/  2n  +  1  /• 


E.  = 


4tt  , _  =J  n  -  1  r; 
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~rr  Pr,  (cos  tt  )  , 


(135) 


I 


where  the  prime  following  the  summation  sign  is  to  indicate  that  the  term  n  =  1 


is  omitted;  the  factor  l/4ff  is  irrelevant, 
give  the  generalized  Stokes'  function: 

Ei  =  ^  [S(rp,  r)  - 


The  series  (135)  may  be  sumn.ed  to 


J_  r  2  _  _5r 

•iff  L  £  r0‘ 


cos  4)  in 


rcosj^  +  i  1 
2r;  J  ’ 


(136) 


see  (Moritz,  1965a, equations  (50)  through  (54)). 

The  density  of  the  generalized  surface  layer  in  this  case  is  denoted  by  X 
instead  of  <p,  so  that  (122)  takes  the  form 


T  =  m  X  E,  dS  .  (137) 

S 

Thenthe  kernel  will  be  given  by  (133a),  multiplied  by  l/4ff,  so  that  (123) 
becomes 


|£  +  S  .JL  rfX  ds. 

drp  r-  4ff  v  rp  r 


(138) 


According  to  (136),  the  main  singularity  of  Ex  is  that  of  1/j l;  the  second, 
logarithmic,  singularity  has  no  effect  in  the  present  situation,  as  will  be  seen 


later.  Hence  the  function  E  behaves  like 

1 

2ff  l 


(139) 


as  l-  0,  that  is,  as  P  approaches  the  surface  S.  Therefore,  the  vertical  derivative 
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of  (137)  will  have  on  S  the  same  discontinuity  as  the  potential  of  a  simple  layer 
(the  function  T  itself  will  be  continuous).  According  to  (110),  this  discontinuity 
is  -  2irocos/3;  in  agreement  with  (139),  X  takes  the  place  of  2 n<p,  so  that  now  the 
discontinuity  is  -  Xcos/3  .  Hence  (138)  becomes  on  S 

fl  +S  CS.  (140) 


Since  the  right-hand  side  of  this  equation  is  equal  to  -Ag,  we  obtain  the  integral 
equation 

Xcos0'S  JJ  A  ds  = 

J 

for  determining  X  from  Ag. 

The  planar  approximation  is 


(141) 


Xcos/3  -  jT X  -  ~ahp  sec ^ da  =  Ag  . 
a 


(142) 


Since  the  kernel  becomes  zero  for  the  sphere  (h  =  h,  =  0),  the  second  term  will  be 


small,  so  that  this  equation  lends  itself  to  an  iterative  solution: 
XG)  =  Agsec/3, 

X^'*  =  sec/3[Ag+~  ;rx(  ^  sec/3daj  . 


(143a) 


(143b) 


We  shall  be  satisfied  here  with  the  linear  approximation 


.  R  -  ,  h  -  hp 

X  '  2?  X~L^  da  ‘ 
o 


of  which  the  solution  is 


(144) 
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1  ■» 


Ag  da  . 


(145) 


A 


h  -  h- 


According  to  (120)  this  is 

X  =  Ag+G,  •  (145') 

Let  us  now  consider  the  linear  approximation  to  Ex  if  P  lies  on  the  surface 
S  itself.  To  this  approximation  we  have  £  =  L  =  2Rsin  (0/2),  so  that  (136) 
reduces  to 


E, 


1^25  3f0 

-  -  cost  -  ^r 


3 

R 


cos  0  in 


R  -  Rcost  +  f; 
2R 


) 


and  finally,  exlcuding  the  zero-degree  harmonic,  to 


E. 


(146) 


where  S(0)  is  the  ordinary  function  of  Stokes.  Hence  equations  (137),  (145'), 
and  (146)  again  lead  to  the  solution  (118)  considered  in  the  preceding  section. 

If  P  lies  outside  the  surface  S,  then  (146)  no  longer  holds  (not  even  as  a 
linear  approximation),  and  (136)  must  be  applied.  In  this  case  we  have 

T  =  —  ff  (Ag  +G1)[s(rP,0,  r)-^—  ]  da  ;  (147) 

a  p 

this  equation  provides  a  simple  and  useful  formula  for  computing  the  external 
gravity  field. 

Hirvonen  (1960)  considered  the  zero-degree  approximation  that  underlies 
this  method:  he  uses  a  function  that  is  essentially  equivalent  to  (136)  but  disregards 
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(Strictly  speaking,  the  kernels  corresponding  to  Ex  and  Eg  contain  a  spherical 
harmonic  of  the  first  degree,  but  this  harmonic  will  integrate  to  zero  in  (138) 
and  (152),  so  that  we  can  disregard  it  here.)  From  the  interpretation  of  (151) 
as  the  potential  of  a  volume  distribution  it  follows  that  neither  T  nor  dT/drP 
will  undergo  a  discontinuity  if  P  crosses  the  surface  S  (p.  57).  Hence  even  for 
P  on  S  equation  (152)  holds,  so  that,  with  the  planar  approximation  R/rP  =  1, 

fj  f  dS  =  Ag  (153) 

S 

is  the  desired  integral  equation  for  determining  n  from  Ag. 

This  integral  equation  is  formally  the  simplest  considered  so  far.  This 
does  not  necessarily  mean  that  it  is  the  most  practical.  Integral  equations  of  the 
first  kind,  such  as  (153),  are  less  tractable  than  those  of  the  second  kind 
(Courant  and  Hilbert,  1953,  p.  159).  (All  of  our  previous  integral  equations 
were  of  the  second  kind.)  The  linear  solution,  however,  is  easily  found. 

As  a  linear  approximation  we  have 
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Comparison  with  (16)  shows  that 


(156) 


so  t^at,  apart  from  a  constant  factor  and  a  small  additional  term,  u  is  essentially 
equivalent  to  the  anomalous  vertical  gradient  of  gravity.  It  is  well  known  that  this 
vertical  gradient  is  much  more  irregular  than  the  gravity  anomaly  itself.  This 
confirms  what  has  been  said  above  about  the  tractability  of  (153),  being  an  integral 
equation  of  the  first  kind. 

It  would  be  possible  to  use  (155)  with  (151),  but  it  is  more  practical  to  con¬ 
sider  the  linear  approximation  of  Ea  rather  than  its  rigorous  expression.  We  have 
r  =  r„  +  ( h  -  hp )  ; 

hence 


Ea  =  & 


dr 


(h  -  h, )  , 


where 


j  _  JL  y>  1 


r  S'  ~  pn  (cos 0)  S'r—  p„  (costf.) 
rP  K-o  n  -  l  K  n  -  i 


corresponds  to  the  strictly  spherical  case.  Differentiating  (148),  setting 
r  =  rP  =  R,  remembering  the  well-known  spherical-harmonic  expansion  of 
Stokes'  function  and  neglecting  very  small  terms  we  find 

=  _L  « /  n  n  •  1  *  2n  + 1  n  S  (ib) 

or  R*  ,1,  —  p-  *  W  ,1.  —  p.  *  • 
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We  substitute 


=  Es0  +2R5”  S<^>  "M  <157> 

in  (151)  and  use  (156),  neglecting  the  second  small  term  which  is  in  agreement 
with  the  planar  approximation.  Thus  we  find 

T  =  R3  JT  n  Ea0  da  +  |  jj’  M  (h  -  hp )  S  (4 )  do 
a  a 

=  R3  n  Ezc  da  -  ~  rp  (h  -  he )  S  (0 )  da  . 
a  a 

The  first  term  on  the  right-hand  side  is  strictly  spherical  and  must  therefore 
be  identical  with  Stokes’  integral: 

R3  JJ  n  Ez°  cfcr  =  ^  jV  Ag  S  (0)  d o  .  (158) 

a  a 

Hence  we  finally  have 

T  S(*,d0  <159> 
a 

as  the  linear  solution  of  (153). 

The  integral  equation  (153),  without  solution,  was  given  by  Brovar  (1964a). 

A  PHYSICAL  INTERPRETATION.  As  Brovar  (1964a) pointed  out,  a  function 
similar  to  Es  can  be  interpreted  as  the  potential  of  a  certain  volume  distribution. 
We  shall  now  give  such  an  interpretation  to  E3  itself  (Brovar’s  E3  is  different 
from  oufs). 
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The  potential  of  a  volume  distribution  is  given  by 

T  ■  f  fj  — ^  -  ’  j  f  ^  r’*  sinQ'de'dA'  dr'  (160) 

r/=o  or 

(the  gravitational  constant  is  taken  to  be  unity).  Assume  the  density  o  to  have  the 
form 

p(r',  e',X/)  =  f(r')MeU').  (161) 

Then  (160)  reduces  to 

T  *  JJ  y  F  da  ,  (162) 

a 


where 


F  =  |  j  f  (r')  r '*  dr'. 


(163N 


Using  (128)  this  becomes 


j  f  (r')  r/B+a  dr'. 


In  order  to  simplify  the  interpretation,  we  assume  that  the  anomalous  potential 
does  not  contain  spherical  harmonics  of  degrees  zero  and  one.  Then  the  sum  in 
(164)  can  start  with  n  =  2.  Furthermore  we  set 


(164) 

T 


f(r') 


.'-4 


Thus  (164)  becomes 


(165) 


F  = 


1 

n  -  1 


P,.  (cos  $  )  . 


Comparison  with  (148)  shows  that,  apart  from  zero-degree  and  first-degree 
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harmonics  left  out  of  consideration. 


c  :  ^  ^ 

r  R  * 


(166) 


Hence  (162)  becomes 


T  *  r  j>  Ead° 

Q 


and,  with  (93b) 


T  =  rr  e2  ds. 
s  R 


(167) 


Comparing  this  with  (151)  we  see  that 
V  =  K4  fi  sec  fi  . 

Hence  the  material  density  of  the  fictitious  volume  distribution  that  produces  T  is 


p  =  psec  fi  ; 

\  I  ✓ 


(168) 


it  increases  towards  the  center. 

Hence  the  function  Ea  may  be  considered  to  correspond  to  the  potential  of 
a  volume  distribution,  and  the  generalized  Stokes’  function  (136)  corresponds 
according  to  (149)  to  the  linear  combination  of  the  potentials  of  such  a  volume  dis¬ 
tribution  and  of  a  surface  layer.  Since  the  potential  of  a  volume  distribution  is 
everywhere  continuous  together  with  its  first  derivatives,  there  is  no  discontinuity 
of  (152)  on  S,  and  the  discontinuity  of  (138)  is  solely  due  to  (139),  which  represents 
the  surface  layer. 
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t  7.  BJERHAMMAR'S  INTEGRAL  EQUATION 

The  approach  of  Molodensky  (sec.  5)  and  its  generalization  by  Brovar  (sec.  6) 
use  surface  layers  on  the  physical  su  rface  of  the  earth.  Bjerhammar  (1964)  attempts 
the  representation  of  the  external  gravity  field  by  a  layer  on  the  reference  ellipsoid 
or,  as  a  spherica*  approximation,  on  a  sphere.  For  this  purpose  it  is  necessary  that 
the  external  gravity  field  can  be  analytically  continued  down  to  sea  level.  This  gives 
rise  to  certain  theoretical  and  computational  difficulties  which,  being  irrelevant  to 
the  linear  approximation,  will  not  be  discussed  here;  see  (Moritz,  1964,  1966). 

The  external  gravity  field  is  thought  to  be  generated  by  a  set  of  fictitious 
gravity  anomalies  Ag*  on  the  sphere  representing  the  reference  ellipsoid.  The 
actual  gravity  anomalies  Ag  on  the  earth’s  surface  as  obtained  by  measurement  are 


then  related  to  Ag*  by  the  usual  "upward  continuation  integral"  (e.g. ,  Heiskanen 


and  Moritz,  1967,  sec.  6-8): 

a 


4g,  =  R*<r>  -  r3)  *£ 


4ffr0 


x-  r 

o 


do  . 


(169) 


The  notations  are  evident  from  Fig.  8.  The  function  Ag  on  the  earth's  surface  being 
given,  this  equation  is  a  linear  integral  equation  of  the  first  kind  for  Ag*.  The 
peculiar  simplicity  of  this  equation  rests  on  the  fact  that  the  integration  is  now 
rigorously  extended  over  a  sphere.  It  was  suggested  by  Bjerhammar  (1964). 

Once  Ag*  has  been  found,  all  computations  can  be  done  directly,  nothing 
but  spherical  formulas  being  involved. 
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Tlie  comparison  of  (173)  and  (172)  provides  an  independent  derivation  of  (16). 

Arnold  (1965)  obtained  the  linear  solution  directly  from  (169)  by  an 
approach  whose  mathematical  justification  is  more  difficult. 

To  get  higher  approximations  for  Ag*,  the  integral  equation  (170)  may 
be  readily  solved  by  an  iterative  method  similar  to  (143a,  b).  As  a  planar 
approximation,  (170)  reduces  to 
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B.  TOE  UNEAR  SOLUTION  AS  A  CONSEQUENCE 
OF  TOE  GRADIENT  SOLUTION 


8.  TOE  VERTICAL  GRADIENT  OF  GRAVITY 
in  fbrt  B  we  stall  derive  the  linear  solutions  obtained  in  Part  A,  and  also 
others,  from  a  almpfe  and  intuitively  evident  formula  which  uses  "free-air 
anomalies  at  sea  level" 


Ag*  =  Ag  -  h  . 

on 

As  a  preliminary  step  we  consider  now  various  expressions  for  the  anomalous 
vertical  gradient  of  gravity,  SAg/Sh.  As  a  spherical  approximation, 
d4g  _  dAg 

dh  dr  ’  (175) 


EXPRESSION  IN  TERMS  OF  Apr, 
which  may  be  simplified  as 


In  sec.  1  we  have  found  the  expression  (16), 


dAg  _ 
dh  2v 


da  , 


(176) 


because  the  small  term  (-2Ag/R),  when  multiplied  by  h  and  subtracted  irom  Ag 
according  to  (173).  will  produce  a  term  of  the  order  of  h/R  and  can  therefore 


be  neglected  as  a  planar  approximation. 


Since  the  integrand  of  (170)  decreases  rapidly  with  increasing  distance, 

we  may  replace  the  sphere  by  its  tangential  plane  and  compute  the  anomalous 
gradient  by  the  integral 


if  *?***-£  i  i-V***- 

a=o  s=o  (177) 

which  is  formally  extended  over  the  whole  plane;  instead  of  x  and  y  also  polar 
coordinates 


s 


=  Vx2  +y2  , 


(178) 


a  =  arc  tan-^-  , 
x 

may  be  used.  The  origin  of  the  coordinate  system  is  at  P. 

EXPRESSION  IN  TERMS  OF  T  OR  C  .  Inserting  (83)  into  (46)  and 

differentiating  with  respect  to  r  gives,  as  a  spherical  approximation, 

d£g.  =  _  daT  _  2  dT  2 

dr  dr*  r  dr  r^ 

This  equation  is  added  to  Laplace's  equation  AT  =  0,  which  in  spherical 
coordinates  r,  <p,  \  takes  the  form 

daT  ,  2  dT  tan <p  dT  ,  1  d8T  .  1  daT  ^ 

d?~  +  r  37’  IT  ^  +  +  3F  =  °' 

The  result,  on  setting  r  =  R,  is 


dAg  _  2T  ,  1  dT  daT  1  daT  \  /17Qv 

dr  R  R  v  dp  d <f>  cos  p  dX  s 

The  first  term  on  the  right-hand  side  can  be  neglected  as  a  planar  approximation 

and  according  to  (35),  with  p  =  90°  -  0,  there  remains 


*  a2T. 


dr 


By  Bruns  ’  theorem  we  have 

T  =  yc  =  gc, 


(180) 


(181) 
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w^C  for  ttae  present  ^rpoee  m.y  fie  identified  ^  ^ 
Hence  (180)  becomes 


(182) 

write  (182)  more  explicitly 


aT  35  G**Z- 

EXPRESSION  in  terms  np^Kfp  We 


as 


^  =  p  (-  tan  p  |-£  +  +  1  aV  \ 

*  R  ^  a?+  ^  -^r )  (183) 

and  express  the  horizontal  derivatives  of  £  m  term*  # 

C  m  terms  of  the  components  £  and 

V  of  the  deflection  of  the  vertical  by  the  relations 


££  =*  -  n*  dr  . 

d<p  .  cr  =  - 


ax 


R?j  cosp  . 


Thus  we  obtain 


=  G(itan0--^  -  - ^3  .  'N 

VR  Rdp  Rcos  <pb\  )' 

R  is  again  convenient  to  introduce  ta%ential  plane  coordinates  by 
Rd<p  =  dx  ,  R cospdX  =  dy  , 


so  that 


(184) 


(185) 


at  is  of  the  order  of  fig.  the  firs,  ,erm  o„  the  fight-hand  side  is 
negligible  as  a  planar  approximation,  and  there 


^S.  = 


remains 


ah 


=  -G 


■i(p  +  »■ 

\dx  d\ 


dy  y  * 


(186) 
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This  formula  was  used  by  Mueller  (1961). 

GENERAL  REFERENCE  SURFACE.  The  vertical  gradient  of  gravity 
itself  is  according  to  a  theorem  of  Bruns  given  by 

|*  =  -  2gJ- 2w3,  (187) 

d/dH  is  the  derivative  along  the  plumb  line,  g  is  gravity,  J  is  the  mean 
curvature  of  the  level  surface  at  the  point  considered,  and  to  is  the  angular 
velocity  of  the  earth's  rotation. 

The  vertical  gradient  of  normal  gravity  is  correspondingly  given  by 

f?  =  -2yj  -  2u2  ,  (188) 


where  d/dh  is  the  derivative  along  the  normal  plumb  line,  and  J  is  the  mean 
curvature  of  the  normal  level  surface. 

We  divide  (187)  by  g  and  (188)  by  y  and  take  into  account  that 


so  that 


Ms  -  is  ay  _ 
aw  "  aw  aw 


-2(J-J> 


Using  (18$)  we  finally  obtain 
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(190) 


'££  =  "2y(J  -  J)  +  20)*^. 

g 

Ike  difference  of  the  mean  curvatures  of  the  actual  and  the  corresponding 
normal  level  surface  can  be  expressed  in  terms  of  f  ,  which  is  the  vertical 
distance  between  these  two  surfaces.  This  is  a  purely  geometrical  problem, 
which  will  be  solved  in  Appendix  II.  "Hie  result  is 

7  -  J  -  -<2J*  -K)  t-iAat.  (191, 

where  K  is  the  Gaussian  curvature. 


Equations  (190)  and  (191)  hold  for  an  arbitrary  reference  surface.  For 
a  sphere,  their  combination  reduces  to  (182),  small  terms  having  been 
neglected. 

EXPRESSION  IN  TERMS  OF  HORIZONTAL  DERIVATIVES  OF  We  may 


write  (176)  as 


'  S  If  do. 


(192) 


since  by  (39),  neglecting  the  second  term  on  the  right-hand  side  as  a  planar 


approximation, 


Ct)  =  t4. 


(193) 


The  right-hand  side  of  (192)  may  be  transformed  by  Green’s  surface  identity 
(34),  so  that  we  obtain 


dAg  R|| 

dh  2 ir 


to  H  D  /  )  do  . 


(194) 


By  (38)  we  have 


D  (tig,  i-  ^  =  -  i 

f°  *.3  *tl> 

because,  as  a  planar  approximation,  R  sin^j 


.  1 

f,  Rd0 


(195) 


=  2  R  sin  =  f0 .  We  further 


,  dAg 

R "  dx  C0Sa  +"fr  sin«. 


(196) 


where  dAg/dx  and  Ag/Sy  are  the  horizontal  derivatives  of  Ag  in  a  north- 
south  and  an  east-west  direction,  respectively,  and  a  is  the  azimuth. 

By  means  of  (195)  and  (196),  equation  (194)  becomes 

3h  “  2tT  JJ  fT*"  (lx*  COSa  +'t^B  sin«)  (197) 

which  expresses  the  vertical  derivative  of  Ag  in  terms  of  its  horizontal  derivatives 
of  first  order. 

Equation  (194)  may  be  further  transformed  by  Green's  identity  (33). 


The  result  is 


=  g-  rr  -£i0Kl 

ah  2»r  JJ  t 


(198) 


which  expresses  the  vertical  derivative  of  Ag  in  terms  of  its  horizontal 
derivatives  of  first  and  second  order. 
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9.  THE  GRADIENT  SOLUTION 


Let  us  consider  the  analytical  continuation  of  the  external  gravity  field  down 
to  sea  level.  The  corresponding  gravity  anomaly  at  sea  level,  which  will  be 
denoted  by  Ag*,  is  related  to  the  surface  free-air  anomaly  Ag  by 

Ag*  =  Ag  -  h  (199) 

phis  hivher  order  terms. 

This  analytical  continuation  was  called  "free-air  reduction  to  sea  level" 
in  (Moritz,  1964,  sec.  6.4). 

Oice  Ag*  is  known,  the  quantities  pertaining  to  the  external  gravity 
field  can  be  computed  by  purely  spherical  formulas.  In  particular,  the 
anomalous  potential  is  given  by  Stokes'  formula  for  external  space: 

T  (r,  0,  X)  =  —  JJ  Ag*  S(r,  0,  R)  do,  (200) 

G 

where  the  notation  of  (136)  is  used. 

This  formula  gives  T  outside  the  physical  surface  S,  and  also  on  S 
itself.  Hence  the  height  anomaly 


is  found  by 

C(r,  °»  X)=il?  JT  Ag*  S(r,  0,  R)  do.  (201) 

o 

It  is  more  convenient  to  use  the  linear  approximation  to  this  formula. 
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Since 


r  =  R  +  h  > 

we  have 

C  (r,  e,  x )  =  c(R.  e,  x>  +!£  h . 

Now,  £  (R,  0,  X)  refers  to  sea  level  and  is  therefore  given  by  the  ordinary  Stokes' 
formula.  Hence,  taking  also  (199)  into  account,  we  obtain 

« h-  <M2> 

o 

This  formula  gives  the  height  anomaly  £  =  £  (r,  0,  X)  at  the  earth’s  surface. 

Since  various  expressions  for  dAg/oh  have  been  given  in  the  preceding 
section,  there  remains  to  determine  d£/dh.  By  differentiating  the  original 
Bruns  formula  £  =  T fy  we  find 


K  =  A  i  bt 

Ch  bh  ^.yJ  y  dh 


so  that  by  (46) 


dh  y  G 


(203) 


Inserting  this  in  (202)  and  multiplying  by  G  we  find 


T  =  5;  <»«> 

O 

which  is  our  basic  formula,  the  gradient  solution. " 

Higher  approximations  to  (200)  may  be  found  by  the  integral  equation 
approach  of  sec.  7,  but  as  we  have  just  seen,  the  linear  solution  (204)  can  be 

i 

derived  in  an  elementary  way. 


REDUCTION  TO  POINT  LEVEL.  The  sea  level  has  no  preferred  position 

in  this  problem;  reduction  to  any  other  level  may  be  used.  If  this  level  has  the 
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elevation  h.  above  sea  level,  then  h  must  be  replaced  by  h  -  h. .  B  we 

P'“  "  '*’•  "*■'**  “  *•  lml  I—*  thro^h  the  gtowid  point  P 

ynder  consideration  ("point  Wei"),  then  h.  .  h,,  a*i  (204,  reduces 

T  *  £  Jlfc  ’  <h-h,>]  S«)d»,  (205) 

etoce  h  in  the  last  tem,  of  (204)  is,  according  to  our  convention,  identical 
with  h„  so  that  h-h,  =  h,  -h,  =0. 

Tbe  advance  (205)  over  (204,  is  its  sintpiicity.  its  disadvantage  is 
that  the  point  level  varies  from  point  to  point. 
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10.  ARNOLD  TYPE  SOLUTIONS 


The  starting  point  is  (205),  the  vertical  gradient  being  expressed  by  (180). 
R 

T  '  4^T  JJ  [Ag-(h-hp)AaT]  S  (0 )  da  .  (206) 

a 


We  write  this  in  the  form 

T  =  4 ^  ff  AgS(0)do+6T, 
a 


(207) 


where 


R  ^ 

6T  =  "iF  ff  'MAa  T  S(0)  da  . 
a 

As  a  planar  approximation  we  have 


(Appendix  I)  and  hence 


(208) 


(209) 


Aa  Tdo  . 


(210) 


By  means  of  Green’s  surface  identity  (33)  this  is  transformed  into 


(211) 


With 
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this  becomes 


6T  *  »  JT  f  5  J,J,<h-h,)D(T.ji-)<to. 

a  ’  °  (212) 

In  agreement  with  (48)  we  have 

D(T,  h)  =*  D  (T,  h)  =  -  y  (£  tan  0X  +77  tan  0a  ) ,  (213) 

and  in  analogy  to  (195)t 


d(t,  iX  i  .  4 

V  *  it  J  c  a<t>  i7 


a  t 

Rd^  • 


(214) 


where 


3T  J  ar  ^ 

RS*  '  r  RS0"  =-V(?  cosa+^suiot) 


(215) 


because  -a£  /Ra$  is  the  radial  deflection  component  corresponding  to  the 
azimuth  a. 

Hence  (212)  becomes 
R3  1 


6  T  =  "  2ff  JT '  Ta  y  ^  ton  h  +rf  tan  &  )  da 

a 


+  ^  rr 


2ff  f. 

a 


7*  y  (£  cos  a  +  77  sin  a)  do  , 


and  using  (209), 


67  =  ~  4n  JT  y  (€  tan  Pi  +  T?  tan  03  )  S  (ip  )  do 
a 


+  JT  7*  y  (£  cos  a  +  r)  sin  a)  dor . 


(216) 
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* 


Inserting  this  into  (207)  we  obtain 


R 

T  4ir  IJ  tAg  ~  y  (£  tan  +  tjtan /3S )]  S(0 )  da 

+  —  (V  h  -  hB 

2rr  ~  £  2  y  (£  cos  a  +  rj  sin  a)  da. 


(217) 


A  slightly  different  form  is  obtained  by  transforming  the  second  member 
of  (212)  using 


(h-h,)  D(J,  -£)  =  5  [(h.MT>_i  j  -TD(ji,  h) 

By  Green's  identity  (34),  using  (193),  we  find 

£  =  jT(h-h»,TAa(-L> 

Cl  '”o  ' 


=  "  lr  JT  T  • 


2 n  i 


(218a) 


Corresponding  to  (214)  we  have 

K  7~  ,  h  )  =  -  -^4 

'  Jo  40  * 


so  that 


"  9*  T D  (  ~  ,  h^do  =  —  (T  Steak  3 h 

2»  .1  ^  ,  ii^dti  ^  Jj  Tdo  _  (218b) 

Hence  the  second  term  in  CUT)  may  be  replaced  by  the  sum  of  (gig., 
and  (218b),  so  that  we  obtain 

T  ~  4tr  J  C  Ag  -  y  d  tan  pl  +  77  tan  )  ]  s  (0 )  d<7 

'  5F  Sj  f?0"h' -"to*  fjj)  Tdcr.  (2I9) 
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I 


For  later  reference  we  note  as  a  by-product  the  relation 


*  (h_lv)  J'i'f?r(h'h^sin^|^)T(to* 


(220) 


which  is  the  deeper  reason  for  the  equivalence  of  (217)  and  (219). 

Equation  (211)  may  be  transformed  in  still  another  way.  It  is  easily 


verified  that  the  relation 


-‘V'D[(I>-MT,  ^]vd(X,  h)-2TD(J-,  h) 

holds.  The  second  term  on  the  right-hand  side  is  transformed  by 

|  JJ  ■  h>  -  -  |  /J i  TAa  hdo  d  .4.  JJ* T Aa  hSW  ,da. 

a  a 

the  first  and  third  term  by  (218a)  and  (218b).  Using  the  expression  so  obtained  for 
6  T  in  (207)  we  find 

R  r\ 

T  =  «  JTWg-  TA,  h)S(tS)do 


2ff  ^  f,3  •2sin<,|i)Tdo. 


(221) 


The  Unear  solutions  (206),  (217),  (219),  and  (221)  are  equivalent.  Three 
of  these,  (206),  (219),  and  (221),  have  been  obtained  earUer  as  solutions  of 
integral  equations:  (107),  (105).  and  (106).  An  inverse  transformation,  leading 

from  the  form  (217)  to  (206),  was  pubUshed  by  Arnold  (1959b).  The  solution 
(219)  was  derived  by  de  Graaff -Hunter  (1960). 
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and  finally 


Gl  =  2»  JJ  Ag  dtr  .  (226) 

a  0 

Hie  combination  of  (223)  and  (226)  constitutes  Molodensky's  linear  solution, 
which  we  have  derived  from  integral  equations  earlier  in  this  report  (sections  5 
and  6).  The  relation  (225)  between  Gx  and  the  vertical  gradient  was  established 
by  Molodensky  et  al.  (1962b)  in  a  somewhat  different  way. 

We  have  presented  this  derivation  here  in  order  to  obtain  Molodensky’s 
solution  as  a  formal  consequence  of  the  gradient  solution  (204).  The  following 
derivation,  however,  is  still  shorter. 

Denote  by  S/  the  level  surface  that  passes  through  the  point  P 
of  the  earth's  surface  S  (Fig.  9),  and  denote  the  anomalous  potential  on  S'  by 
T#  and  on  S  by  T.  Then,  for  two  corresponding  points  on  the  same  vertical 
such  as  B  and  B*  we  have 

T  =  T'+Th  <h‘h'> 


*  T'+(-Ag-|f)<h-h,) 


or,  as  a  planar  approximation, 

T  =  T'  -  Ag  (h  -  h, )  . 

This  is  substituted  into  (224)  with  the  result 


Ag  +Gi 


n3  '“p f  _  «X*/  n2  u  _  u 

*  -x  - 


2ir 


fc 


2tr 


T"  Ag  do  . 
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I 


i 

I 

i  i 

i 

i  I 

! 

The  first  two  terms  on  the  right-hand  side  give  the  gravity  anomaly 

|  ; 

Ag'  on  the  level  surface  S'.  Since  the  point  P  under  consideration  lies 

!  \ 

both  on  S  and  on  S',  we  have  for  this  point  Ag'  =  Ag.  Hence  the  first 
two  terms  on  the  right-hand  side  will  cancel  with  the  first  term  on  the  left- 
hand  side,  and  there  remains  the  relation  (226)  which  was  to  be  derived. 


12.  TRANSFORMATIONS  OF  MOLODENSKY'S  SOLUTION 
By  means  of  (193)  we  may  write  the  Molodensky  correction  (226)  in  the  form 

Gi  =  Iff  JJ*  ^  ■ hp)  ^As  (r)*  *  <227) 

a  0 

which  lends  itself  to  transformation  by  means  of  Green's  identity  (34).  The 
result  is 

°i  =  "If  jT  D[(h -h„)Ag,  ^-]do. 

In  agreement  with  (195)  this  is 

=  f?  JJ  ^  [<b  “  *V  >  j  &  •  (228) 

This  is  an  alternative  form  of  (226).  Its  interest  is  mainly  theoretical, 
because  it  serves  as  a  starting  point  for  further  transformations. 

By  differentiating  the  product  between  brackets  we  find 

Gi  s  G1X  +Gja  ,  (229) 

where 


The  term  Gia ,  when  inserted  in  Stokes'  integral,  gives  rise  to  a  term 

T>=  "5  J/Oi.  f  Si  ^  “a,  (232) 

a  a  0 


according  to  (209). 
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Being  satisfied  with  the  planar  approximation,  we  may  in  these  small 
correction  terms  extend  the  integration  formally  over  the  infinite  plane. 

Introducing  plane  polar  coordinates  (178)  and  using  the  notations  of  Fig.  10  we 
have 

2ir  00  2*  00 

Tia<0)  J  J  ^a^-!  s  dsdot=~  J  J*  Gia(s,  or)ds  dot , 

0=0  6=0  “=0  85,0  (233) 

,  2ff  •  . 

Gi*(s,  a)=—  J  J  jr  57'  Ag(s'»  or')  s'ds'dor'.  (234) 


a  =o  sf=o 


Here 


consequently 


£08  =  s2  +  s'2  -  2s  s'  cos  ( or '  -  or) ; 


9h  _  3£,  3h  1  9L,  3h 
3£0  9s  9s'  s 3a:  dor' 


(2J5) 


becomes 


3h  _  s/-sco8(a '  -  a)  3h  s  sin  (a *  -  or )  9h 

L  3s'  £„  s'3a'  ’ 


(236) 


Substituting  (234)  into  (233)  and  taking  (236)  into  account  we  obtain 

T's(0)  =  W  .r  !'  /  J  j4{[s'-scos(a'-a)]  + 

/wsn  csrt  /V  o  srt  9 


or=o  s=o  a  '=o  s  =o 


+  ssin(a/  -  or)  77^ j  Ag(s',  a')  s'ds'dc/  •  ds  da  .  (237) 

We  shall  now  perform  the  integration  over  s  and  or  first;  with  respect  to  this 
integration,  Ag,  9h/3s',  and  3h/s'3c/  '.re  to  be  considered  as  constants.  By 
standard  methods  of  integration  we  find 
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'Sr  - :  r - V 

*  *  r  ’  »  ■ 


i 


r  ds 

J  V" 

8=C  # 


00 

=  J*  [s*  +  s,a  -  2s  s '  cos  (a '  -  a)  T3/s  ds  = 


1 


«  “a 
is  sin  — - — 


and 


1  ¥*- 

8=^0  0 


1 


so  that 


qg 

clr/  /»  v  ■»  .  11-  cos  (a '  -  a )  1 

7TT  C*  “  8  '-'0R  -  a)  J  ds  =  — t  - - L  *=  — /  , 

i  l0  8  n  .  a  a  -a  s' 


2  sin 


2ir  ®  z  /  /  v 

n  f  s'-tcos  a  -a  .  .  2tr 
J  J  - TT - L  dsda  =  -r  . 


a=o  s=o 


(238) 


Similarly  we  find 


J  J  ‘"*ft ':R>  **  =  o 


a=o  s=o 


(239) 


by  performing  the  integration  with  respect  to  a  first. 

In  view  of  (238)  and  (239),  equation  (237)  reduces  to 

Tia(0)  =  jjj  \  $  jr  6g(s',ct')s' ds'  da'  . 

a  r=o  s' =o 

If  we  now  return  to  the  sphere,  s  ds'  da'  becomes  R3da,  and  s'  becomes  £0, 
so  that  we  obtain 

Tls  "  2ff  £7  *8  O?  da  JT  ^8  S(^)do,  (240) 


according  to  (209), 
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Substituting  (229)  into  (223)  and  taking  (240)  into  account  we  find 

T  HQ*  +G'i  +  **  r|^)  S(Vp)  dor 

a  y 

or,  to  the  same  accuracy, 

T  -5  JJ'W8+Ci>>0+r|t)s<*)‘‘»-  (241) 

(7  * 

The  term  dh/R3$  is  the  radial  inclination  of  the  terrain. 

An  expression  for  Gn  has  been  found  in  (230).  It  may  be  transformed 
by  means  of  (220),  in  which  T  is  replaced  by  Ag.  We  readily  obtain 

Gl1  =  2ff  (h  *  h>  “sin^ff)  *  (242) 

a  0  * 

Hence,  (241)  with  Gx  expressed  by  (23G)or  (242)  constitutes  another 
linear  solution  of  Molodensky's  problem.  A  form  essentially  equivalent  to  the 
combination  of  (241)  and  (242)  was  obtained  by  Brovar  (1963a,  1964b)  by  solving 
an  integral  equation. 


13.  A  CONNECTION  WITH  TOE  TERRAIN  CORRECTION 
Pellinen  (1964)  suggested  a  solution  of  the  form 

T  JJ  (Ag+G')S«,)do  +  t. 


(243) 


where 


O'  -S- 


(244) 


and  t  is  a  correction  term,  which  will  be  considered  in  what  follows. 

It  will  be  useful  to  denote  by  the  subscript  O  the  point  at  which  T  is  to  be 
computed,  and  to  use  the  subscripts  1  and  2  for  distinguishing  the  variables 
of  integration.  Hence  Molodensky's  formuias  (223)  and  (226)  may  be  written 


Ta  =  ~  JJ  CAgx  +  (G3.K  3  S(^0l)dq  , 


<Gx)i 


~  f*  l  ^8-  ~  Aga  d Os  • 

2n  J-j  fX23  ^  2 

Os 


(245) 


(246) 


The  advantage  of  this  new  notation  is  that  now  (246)  can  be  substituted  into  (245) 
without  danger  of  confusion: 

To  =  ~  Jj(Agx  +“  JJ  Afe  dq8)s(^oX)dcy1  .  (247) 

0X  0a 

It  is  readily  verified  that  the  following  equation,  given  by  Pellinen  (1964), 
is  equivalent  to  (247): 

e>2 

01  03 

•  TSF  JJ-  J’J’ **'  «Oi  <**  .  (248) 

*13 


T”  ■  £  /![«»  +|  ]  sc*.,  )da, 

n,  ru  12 


01  03 
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It  is  of  the  form  (243).  Unfortunately  the  expression  for  t  --  the  last 
integral  on  the  right-hand  side  --  is  impracticable  since  S($0l)  and  S(^l2) 
are  singular  at  the  point  G,  and  there  is  apparently  no  direct  way  of  avoiding  or 
neutralizing  this  singularl  y. 

Therefore  we  shall  proceed  differently.  We  consider  the  gradient 
solution  (204),  the  vertical  gradient  being  expressed  by  (176).  Using  our 
present  notation  we  may  write  it 

To  =  Tff  J^gi  *  |f  J7hi  dffs  )  S  (ill  2 !)  dff!  -h0  Ag0  . 

ffi  CTa  12  (249) 

We  form  the  arithmetic  mean  of  the  equivalent  equations  (247)  and  (249)  and 

perform  some  elementary  manipulations.  The  result  is 


<Tl 


Os 


+  ^r(Agl  N  dffa  '  hl  n  ^  )1  s (0 

Os  13  ff2 


Os 

-  |  h.  Age  . 


(250) 


Tliis  equation  has  the  form  (243),  the  correction  term  t  being  given  by 
1  =  &  JK  S(0)  dCT  -  ih  Ag  ,  (251) 


where 


G"  =  *8  |f  JJ  cla  "  h4|r  SI  >  (252) 


o  -  a 

these  equations  are  written  in  our  old  notation. 
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Although  the  new  expression  for  t  does  not  contain  singularities,  it 
does  not  look  vexy  promising;  still,  under  certain  conditions  of  practical 
significance  it  reduces  to  a  surprisingly  simple  form. 

AN  ADDITIONAL  ASSUMPTION.  Experience  shows  that  the  free-air 
anomalies  Ag  can  often  be  represented  by  a  linear  relation 

Ag  -  a  +  bh  (253) 

where  a  and  b  are  approximately  constants.  In  addition,  b  is  often 
approximately  equal  to  the  Bouguer  gradient 

b  =  2irkp  -  0. 11  mgal/meter  (254) 

(k  =  gravitational  constant,  p  =  density),  so  that  a  is  essentially  nothing  else 
than  the  Bouguer  anomaly,  which  in  this  case  is  largely  independent  on  local 
irregularities  of  topography.  In  statistical  terms  this  expresses  a  correlation 
of  the  free-air  anomaly  with  elevation. 

By  substituting  (253),  with  constant  a  and  b,  equation  (252)  reduces  to 


h  -  hp 


G  ■  4?  a  I J  TT-  . 


(255) 


so  that 


i  G'  ‘  5  JJ  17*  * 


To  the  same  accuracy  we  have  as  a  planar  approximation 

c”  -  -  —  -  —  rr  ^  ~ Jk  a, T 

: G  r  s  JJ  tx  47  • 
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By  (5)  and  (11),  the  inversion  of  this  equation  is  Stokes'  formula: 

h  -  •  \  -  if  If  O' «(*)<*>. 

a 

so  that 

5  IJc'SOOJdo-  =  -§  h. 

a 

Since  Stokes'  integral  is  known  to  suppress  the  zero-degree  spherical  harmonic, 

this  equation  will  hold  only  if  h  does  not  contain  such  a  harmonic,  that  is  if 

J'jhdo  =  O. 
a 

If  this  is  not  true,  then  we  must  obviously  subtract  the  mean  elevation 

"■  =  JJ h  * 

a 

from  h,  so  that  r  finally  have 

£  J/G'  S«,)dj  =  -§(h-M  .  (256) 

a 

Hence  (251'  becomes 

t  =  -  |  a  (h  -  U  )  -  |  h  Ag  t  (257) 

Usually  the  zero-degree  harmonic  of  Ag  is  assumed  to  vaMsh.  By  (253),  this 
gives  the  condition 

°  =  4?  =  a  +bh.  , 
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from  which 

a 

and 


-bh. 


(258) 


Ag  =  b  (h  -  h, ) .  (259) 

Then  (257)  reduces  to  the  very  simple  form 

t  =  -  |  Ag  <h  -  h. ) ,  (260) 

so  that  (243)  becomes 

T=~  JJ  (Ag+G')S(<j>)d0--!  Ag(h-tO-  (26D 

a 

Let  us  now  consider  the  quantity  G '  under  the  assumption  (253). 
Substituting  this  relation  into  (244)  and  expressing  b  by  (254)  yields 

G'  =  \  kp  Ra  JT  lh  do  ,  (262) 

O  0 

which  may  be  shown  to  be  essentially  identical  with  the  conventional  terrain 
correction  for  "deviation  from  the  Bouguer  plate."  This  interpretation,  which 
was  given  by  Pellinen  (1964),  furnishes  an  important  link  with  conventional 
methods. 

Obviously  the  quantity  a  in  (253),  being  the  Bouguer  anomaly,  is  not 
a  true  constant  for  the  whole  earth,  as  was  assumed  in  the  derivation  of  (260), 
However,  for  (256)  to  hold  as  a  planar  approximation  it  is  sufficient  for  a  to 
be,  loosely  speaking,  a  local  constant,  or  rather  a  quantity  that  varies  relatively 
slowly.  Then  h,  is  riot  the  average  height  over  the  whole  earth,  but  rather  some 


88 


’Z* 


V 


» 


local  average,  which  we  may  consider  to  be  defined  by  (258)  where  a  is  the 
Bouguer  anomaly  and  b  is  the  Bouguer  gradient. 

To  get  an  estimate  of  t,  consider  its  effect  on  the  height  anomaly 
£  =  T/G.  This  effect  is  by  (260) 

*  G1 

or  by  (259) 

6C  =  35  (h  -  hj2.  (263) 

If  h  -  =  1000  km,  which  is  a  very  extreme  case,  then 

=  5cm. 

Hence,  in  those  cases  where  Ag  satisfies  a  relation  (253)  to  a  sufficient 
accuracy  and  extent  for  (261)  to  hold,  we  can  probably  always  neglect  the  term 
t  altogether  and  use  the  formula 

If  WS+G')  S(«do,  (264) 

a 

that  is,  adding  to  Ag  the  terrain  correction  G'  instead  of  the  Molodensky 
correction  Gx . 

This  solution  is  of  particular  interest  from  the  point  of  view  of  deflections 
of  the  vertical;  see  sec.  15. 
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C.  DISCUSSION  AND  APPLICATIONS 


14.  HEIGHT  ANOMALIES 

In  the  preceding  Parts  A  and  B  we  have  found  expressions  for  the  anomalous 
potential  T  at  the  physical  surface  of  the  earth.  These  formulas  give  directly 
the  height  anomaly  £  as  well,  since  by  Bruns'  theorem 

£  =  f  (265a) 

or,  as  a  spherical  approximation, 

£-§.  (265b) 

where  y  is  normal  gravity  and  G  an  average  value  of  y  for  the  whole  earth. 

We  shall  now  collect  the  main  expressions  for  T  previously  found,  it 
being  understood  that  by  dividing  them  by  the  constant  G  according  to  (265b) 
we  obtain  expressions  for  £  . 

A.  GRADIENT  SOLUTIONS  (sections  8  and  9) 

*•  t  =  ^  h!) 
a 

«•  T  *-£  J'J’fc- Th  <"-M]s«>da 

a 

Vertical  gradient 

a)  by  measurement 

b>  iaa  =  £  rr  \Asr  * 

'  ^h  2lt  JJ  £„3 
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*>**•«! 


c)  =  A  T  _ 
3  h  2 


*>  *g) 

8  h  2ir  ff  J75"  C  cos  “  +  sin  a )  da 

«  Ih  =f  Jj‘^  da 
a  ° 

B.  MOLODENSKY  TYPE  SOLUTIONS  (sec.  11) 

IJ 

UI-  r  =  Tv  JJ  (^g  +G i  )  S  (0)  da 

a 

a>  °>  JT^  Ag* 

a  *9 

b>  Cx  =  ~  rr  3!£jL  if.  ■  v  1 

2n  JJ  £o3  L^1  h»^^gJ  dcr 

C.  BROVAR  TYPE  SOLUTIONS  (sec.  12) 

IV'T  =  4a  JJ<‘3S+Gi.)(l+^)s«,)da 

a>  G»  =|  //$**-*>$  da 

b)  G'“  '  2a  JJ  «7\h  ‘  h>  '  sin  )  ag  da 
D.  ARNOLD  TYPE  SOLUTIONS  (sec.  10) 

V*  T  ~  4?  JT^-yUtan/^  +T?ton/Sa)]S(0)da  +  K 


ca2c 


Without  exception,  these  solutions  consist  of  the  original  Stokes' 
integral  and  a  correction,  which  assumes  very  different  forms.  We  shall  now 
briefly  discuss  these  forms. 
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As  for  the  gradient  solutions,  the  basic  difference  between  I  and  II  is  that 
I  uses  sea  level  (or  any  other  fixed  reference  level),  whereas  II  uses  a  point 
level  that  varies  from  point  to  point.  Form  II  is  somewhat  simpler,  but  I  is 
more  suitable  for  large  scale  data  processing. 

The  vertical  gradient  dAg/Sh  entering  in  these  formulas  may  be  obtained 
from  measurements  (a),  but  the  respective  techniques  and  instruments  are  still 
at  the  experimental  stage.  The  most  practical  form  of  computing  it  is  probably 
from  gravity  anomalies  using  formula  (b).  Equation  (c)  expresses  it  in  terms 
of  the  disturbing  potential  T  or  the  height  anomaly  (or,  approximately,  die 
geoidal  height)  £,  but  since  second  horizontal  derivatives  are  involved,  the 
knowledge  of  £  must  be  so  detailed  and  accurate  that  the  use  of  this  equation 
is  hardly  feasible  in  practice.  Equation  (d),  which  expresses  dAg/dh  in  terms 
of  the  components  £  and  r\  of  the  deflection  of  the  vertical,  is  more  useful 
practically;  even  astro -geodetic  deflections  may  be  used.  Form  (e)  requires 
as  data  the  horizontal  derivatives  of  Ag,  which  could  in  principle  be  directly 
obtained  by  torsion  balance  measurements.  Otherwise  in  (e),  just  as  in  (f), 
the  gravity  anomaly  Ag  could  be  used,  but  these  two  expressions  are  definitely  less 
practical  than  (b). 

Molodensky’s  solution  III  is  very  simple  and  practical  when  Gj  is 
expressed  in  the  form  (a);  (b)  has  only  theoretical  interest. 

An  advantage  of  the  solutions  of  type  IV  over  those  of  type  HI  is  that 
the  integrals  entering  in  Gu  are  better  convergent  than  those  entering  in  Gx . 
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Mathematically  speaking,  the  integrals  IV  a,  b  are  absolutely  convergent,  whereas 
die  integrals  III  a,  b  only  converge  conditionally  (Bxovar,  1963a).  Ibis 
advantage  is  made  up  by  the  fact  that  IV  is  more  complicated  than  HI,  the 
inclination  of  the  terrain  entering  directly. 

The  Arnold  type  expressions  have  a  certain  theoretical  significance, 
since  they  are  essentially  solutions  of  the  original  integral  equation  for  the 
anomalous  potential  T  (sections  2  and  4).  The  forms  V  a  and  VI  have  a  certain 
esthetic  priority  since  they  express  the  correction  to  Stokes'  integral  only  in 
terms  of  £  and  tj  or  only  in  terms  of  T  ,  respectively.  The  practical  signifi¬ 
cance  of  these  expressions  is  small.  (We  have  chosen  the  name  "Arnold  type 
solutions”  because  Arnold  (1959a)  was  among  the  first  to  find  one  of  these 
solutions;  later  (1959b)  he  transformed  this  into  the  gradient  solution,  which 
he  recommends  for  practical  use.) 

The  interest  of  the  solutions  of  Pellinen's  type  rests  in  the  fact  that  VII 
contains  the  expression  G '  which  (1)  converges  better  than  Molodensky's  Gx 
and  (2)  has  a  close  relationship  to  the  conventional  terrain  correction.  How¬ 
ever,  the  general  expression  (a)  for  the  correction  t  is  impractical.  Only  in 
the  case  that  the  gravity  anomalies  are  strongly  correlated  with  elevation  does 
this  expression  reduce  to  the  simple  form  (b);  in  this  case,  furthermore,  the 
term  t  may  even  be  completely  neglected  in  practice.  Although  this  condition 
for  (b)  to  hold  will  very  often  be  fulfilled  in  practice,  it  should  be  kept  in  mind 
that  VII  b  is  not  universally  valid  to  a  linear  approximation  as  all  the  other 
solutions  are. 
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As  a  conclusion,  for  practical  application  we  are  left  with  die  gradient 
solutions  1  and  II,  the  gradient  being  computed  by  (b)  (in  exceptional  cases, 
the  use  of  (d)  might  possibly  be  feasible)  and  with  Molodensky's  solution  III  a. 
In  spite  of  formal  similarity,  the  solution  II  is  not  of  the  form  111,  since  G  x 
is  a  functior  only  of  its  position  on  the  earth's  surface,  whereas 

-  <h-h,) 

depends,  in  addition,  also  on  the  computation  point.  For  the  same  reason,  11 
is  somewhat  inferior  to  1  and  111  with  respect  to  large-scale  computations; 
however,  this  form  is  eminently  suitable  for  deflections  of  the  vertical,  as  we 
shall  see  in  the  next  section.  A  drawback  of  III  a  is  the  strong  and  direct 
dependance  of  Gi  on  elevation,  in  particular  on  inclination  since  (h  -  hp  )/f0 
is  of  the  order  of  inclination,  Hus  drawback  is  avoided  in  the  slightly  longer 
form  1,  since  1  b  does  not  contain  the  elevation  as  III  a  does. 

Hie  Brovar  type  solutions  IV  a,  b  show  certain  interesting  features 
which  make  them  worth  to  be  further  investigated,  although  they  seem 
definitely  less  practical  than  Ilia. 

The  solution  VU  deserves  further  study  in  view  oi  its  relation  to 
conventional  methods;  in  the  following  section  we  shall  see  that  it  has  also 
particular  interest  for  computation  of  deflections  of  the  vertical. 
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15.  DEFLECTIONS  OF  THE  VERTICAL 
According  to  equation  (16)  of  (Moritz,  1964)  the  components  £,  tj  of  the 
deflection  of  the  vertical  are  given  by 


=  .  K  .  & 

dx  y 


tan  /31  , 


(266) 


n  =  ’  ^  tan  » 


where,  as  usually , 


a  dh  .  dh 

tanft  =  —  ,  tanfc,  =  —  . 


(267) 


the  derivatives  b/bx.  and  b/by  being  taken  along  the  local  horizon  in  a  northern 
and  an  eastern  direction. 

Each  deflection  component  thus  consists  of  two  terms:  (1)  a  derivative 
of  the  height  anomaly  and  (2)  a  term  depending  on  the  inclination  of  the  terrain. 

The  derivatives  b£/bx  and  /by  are  found  by  differentiating  any  of  the  solutions 
listed  in  the  preceding  section.  We  shall  here  limit  ourselves  to  the  most  important 
types  A,  B,  and  E;  the  other  types  may  be  treated  in  the  same  way. 

Let  us  start  with  Molodensky's  solution  III: 


J7  (Ag  +  G:  )S  (ip)du. 


(268) 


Differentiation  of  Stokes'  integral  gives  Vening  Meinesz'  formula,  to  which 
according  to  (266)  the  inclination  terms  must  be  added.  The  result  is 


- - —  V.VL  S 


in  the  second  tern.,  y  has  been  replaced  by  n.ean  gravity  G.  An  equivalent 
formula  is  found  in  (Molodenskii  et  al. ,  1962a,  p.  124). 


In  the  second  place  we  consider  the  gradient  solution  I: 
t  =  4^G  S(*)d,-G  hi8' 


(270) 


Differentiation  of  the  second  term  on  the  right-hand  side  of  this  equation  and 
use  of  (267)  give 

1  .  .  \  _  i  ,  ^Ag  1  ^h 

Co  hAs )  ‘  "  h  +  ~ 


G  ^  \  G  ^  ^ 


Jl  4.  a 

-  G  77  +  G  ,an  01 


Hence  we  obtain 


5  *  JTg  j"  Q*  "fh  h)  7*  cos  "do  +  g  f  0.  • 


The  inclination  term  cancels  out,  and  there  remains 


{{}  ■  470  Q  *+H  &&}•«»> 


Let  us  now  compare  (269)  with  (271).  Theoretically,  these  two  solutions 
are  completely  equivalent.  This  does  not  mean,  however,  that  they  are  equally 
suited  for  practical  application.  In  both  cases  the  correction  terms  may  assume 
large  values.  As  an  example,  consider  (269),  assuming  Ag  =  100  mgal  and 
0i  -  45°.  Then 


tan  0!  *  1CT4  =  20"  . 


Tb  have  a  comparable  situation  for  (271),  consider  an  elevation  h  =  1000  m  and, 
again,  an  inclination  =  45°,  assuming  a  linear  dependence  of  the  free-air 
anomaly  on  elevation, 

Ag  =  a  +b  h  , 

with  b  =  0. 1  mgal/m.  Then 

dAg  ,  dh  . 

h-  =  b,ar.ft  , 

so  that  in  (27 ! ) 

^  =  10-4  =  20* 

G  dx 


as  before. 

The  essential  difference  between  (269)  and  (271)  is  that  the  correction 
term  in  (269)  is  independent  of  elevation,  whereas  the  correction  term  in  (27 1) 
is  proportional  to  elevation,  so  that  by  a  suitable  choice  of  the  reference  level 
it  can  be  made  as  small  as  we  like.  By  reducing  the  free-air  anomalies  to 
point  level  (sec.  9)  it  becomes  zero;  then  (271)  becomes  simply 


This  result  could  also  have  been  obtained  by  formal  differentiation  of  the 
gradient  solution  II,  but  the  mathematical  justification  of  this  procedure  is 
not  immediately  obvious. 

Summarizing  we  may  say  that  the  correction  term  in  both  (269)  and  (271) 
show  an  undesirable  dependence  on  the  inclination  of  the  terrain.  In  (269),  this 
dependence  is  direct;  in  (271)  it  is  indirect,  through  an  approximate  linear 
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dependence  of  the  free-air  anomalies  on  elevation.  In  (269),  the  effect  of 
inclination  cannot  be  removed  or  mitigated  in  any  way;  whereas  in  (271)  it 
can  be  made  smaller  by  reckoning  the  elevation  from  an  average  level 
approximating  the  terrain  rather  than  from  sea  level,  and  it  can  be  made 
zero  by  reckoning  the  elevation  from  point  level,  thus  arriving  at  (272). 

I 

The  correction  to  the  simple  Veiling  Meinesz’  formula  for  £  is 
according  to  (272)  expressed  by 


4ffG"  H  fh*  <h'h')  C0S“da  • 
a 


d  0 


(273a) 


according  to  (269)  by 


Jr.  IS 


4flG 


d  S 
d  ^ 


cos  a  da 


tan  Pi  • 


(273b) 


and  according  to  (271)  by 


1  pp  d  S  .  ,  h  dAg 

—  I  I  -f*  h  -7  coso;d(T+  -  -7*  . 
t  G  ^  h  di/i  Gdx 


4ffG  JJ  3h  d  ib 
a  Y 


(273c) 


This  correction  is  computed  by  (273a)  as  a  single  small  term,  whereas  in  (273b) 
and  (273c)  it  is  obtained  as  the  difference  of  two  larger  terms  computed  in  a 
different  way  and  consequently  affected  by  errors  which  may  seriously  endanger 
the  result. 

TRANSFORMATION  OF  MOLQDENSKY'S  FORMULA.  It  is,  however, 
possible  to  transform  Molodensky's  formula  (269)  in  such  a  way  as  to  eliminate 
the  computational  problems  caused  by  the  explicit  occurrence  of  the  inclination. 
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whence 


or,  by  (176), 


Using  the  notation  of  sec.  13,  we  equate  (273a)  and  (273c): 
4ffG  jfJ*  cosacbi  = 

s  ’  *5  {{ (!£) hi  if  “s«d®+ Kf£)0 . 

G  (|? ).  =  4ffG  /i*  ^  (jh)x  5f*  cosadtTi 


^  9  4irG  J‘J'L2ff  W  KI7  '^)dcral  j^cosado, 

ffj  %  B  J  Oip  1 


ms  identity  holds  for  arbitrary  functions  h  and  tig.  We  may.  therefore, 
interchange  h  and  tig;  wtth  ah/Sx  =  tan  ft  we  this  obtain 

=jjg  /;[£  jt  £  Ag.  h  -h,)dohi  «  cos addr 

ffi  Oa  12  -IflV 

TTus  expression  is  substituted  in  f273hl  anH  n  i.  . , 

n  V/i3D)>  and  G  is  expressed  by  (246).  The 

result  is 

«5  Jj" [lr  Jj  <d&  -  Ag„ )  (It  -  hjdts  ]  —  costrdo  , 

which  is  a  correction  term  to  be  added  to  the  simple  Vening  Meinesa  formula 
for  £,  it  is  equivalent  to  any  of  the  expressions  (273). 


Hence  we  obtain,  in  our  usual  notation, 

0 =  «c  {^s“j  do, 


(274) 


where 


g  ,  =  r  r 11 " 


*>■  pp  n  - 

2ir  J  J  “  o  3  (Ag  -  £gA )  c b 

*o 


(275) 
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is  a  new  correction  term  which  differs  from  the  Molodensky  correction  G  ^  by 
Ag  -  Ag»  taking  the  place  of  Ag.  Here  AgA  is  the  gravity  anomaly  at  the  point 
at  which  £  and  r)  are  to  be  computed,  whereas  Agp  would  be  the  anomaly  at 
the  point  to  which  G,  refers. 

USE  OF  TERRAIN  CORRECTION.  We  shall  finally  consider  the  solution 
VII  b  of  the  preceding  section: 

C  =  ~  J/(Ag+G')  5  (0 )  dcr  -  2jjr  Ag  (h  -  h. ) ,  (276) 

a 

which  holds  when  there  is  a  sufficiently  strong  correlation  of  the  free-air 
anomalies  with  elevation;  then 

c>  _  pp  (h  -  hp)  (Ag  -  AgP) 

4tr  {3 

n  a 


This  term  cancels  precisely  with  the  second  term  on  the  right-hand  side  of  the 
first  equation  of  (266),  so  uutt  Uiere  remains 
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(277) 


nJ  4fG  JJ  “  7  d  A 

o 

REVIEW.  We  shall  now  collect  the  results  so  far  obtained  and  classify 
them  in  the  same  way  as  die  solutions  of  sec.  14,  but  prefixing  the  letter  D  to 
the  roman  numerals,  so  that  the  solution  D1  for  deflections  corresponds  to 
the  solution  I  for  T  or  £  . 

A.  GRADIENT  SOLUTIONS 

“•  < =  «g  |r(4«‘l^h) H cosodff+! 


DU-  « =  So  JJfc'ff  (h'h>,]di c°8ad0 


diu. 


B.  MOLODENSKY  TYPE  SOLUTIONS 


•  «  =  J/<a8+Gi)  df  cos“d<’-lf 


d  4 > 


where  G,  =  jj-  JJ  da 


f  =4?G  JJ' >8 +G>  iff  cos  a  da 


d  0 


where  G,  =  f  ft  ~  h».>  <*  ~  da 
1  2ir  p. 3 


DVII 


E.  PELLINEN  TYPE  SOLUTION 


b*  *  =  4?G  JT  (Ag+G,)^f  cos  ado 


where  G'  =  —  jjO1  -  M  (Ag  -  AgP )  d(y 
4tt  ^  £s3 


102 


As  compared  to  Sec.  14,  the  present  collection  contains  not  all  solutions  given 
there  but  only  those  which  are  believed  to  be  most  significant  for  the  computation  of 
deflections*  Only  expressions  of  the  deflection  component  £  are  given,  the  corre¬ 
sponding  expressions  for  rj  being  evident.  Here  A  is  the  point  at  which  £  and  r\  are 
to  be  computed. 

For  practical  application  those  solutions  that  do  not  contain  the  inclination 
tan  /3:  or  the  horizontal  derivative  oAg/Sx  are  preferable  for  computational  reasons 
given  above  (a  small  term  should  not  be  computed  as  the  difference  of  two  larger 
values  obtained  in  different  ways).  Therefore,  DII  (reduction  to  point  level)  is 
preferable  to  D  I,  and  D  I1I#  is  preferable  to  D  III. 

A  comparison  of  the  expression  (179)  for  dAg/bh  with  that  of  Gx  given  above 
shows  readily  that  the  gradient  solution  D  11  is  computationally  simpler  than  D  III'. 

In  addition,  Gx  and  Gj  depend  strongly  on  the  irregularities  of  topography. 

The  Pellinen  type  solution  D  VHb,  when  applicable,  shares  the  computational 
advantages  of  D  11  and  D  III7.  Moreover,  G '  is  a  true  function  of  position  just  as  Gx , 
having  a  unique  value  at  every  point  of  the  earth's  surface  (and  being  consequently 
representable  by  a  map  or,  say,  by  a  set  5'x  5'  mean  values),  whereas  (dAg/dh)(h-h*) 
and  depend,  in  addition,  on  the  computation  point  A.  On  the  other  hand,  the 
applicability  of  D  I,  D  II,  D  III,  and  D  HI'  is  unrestricted  (as  tong  as  the  linear 
approximation  is  sufficient),  whereas  D  Vlib  presupposes  a  strong  correlation  of  the 
free-air  anomalies  with  elevation. 

Formulas  of  a  different  type  will  be  given  in  sec.  17. 


16.  SPHERICAL  HARMONICS 


The  conventional  formulas  for  computing  coefficients  of  spherical 
harmonics  of  the  earth's  external  gravitational  potential  require  the  free-air 
gravity  anomalies  to  be  given  on  a  sphere.  As  far  as  the  anomalous  gravity 
field  is  concerned,  it  is  permissible  to  identify  this  sphere  with  the  reference 
ellipsoid  —  this  is  the  snkerical  approximation  —  but  not  with  the  physical 
surface  of  the  earth  to  which  the  free-air  anomalies  primarily  refer,  because 
the  inclination  of  the  terrain  is  not  negligible;  cf.  sec.  3. 

The  free-air  anomalies  at  sea  level  as  defined  in  sec.  9, 

Ag*  -  Ag  -  h  , 

however,  are  directly  suited  for  computing  the  spherical -harmonic  coefficients 
by  means  of  the  spherical  formulas,  since  they  refer  to  the  reference  ellipsoid 
which  is  represented  by  a  sphere. 

Let  Ag*  be  a  Laplace  surface  harmonic  so  computed  from  Ag*,  then 
the  harmonic  T„  of  the  same  degree  n  of  the  a^malous  potential  will  be  given 
by 

Tn  =  Ag*  .  (278) 

Hence  the  gradient  solution  (free-air  reduction  to  sea  level)  is  the  most  direct 
method  for  computing  spherical  harmonics. 

Let  us  now  consider  Molodensky's  solution  (sec.  11).  The  connection 
between  the  anomalous  gradient  and  the  Molodensky  correction  G;  is  according 
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to  (225)  given  by 


.  _  r  Rf.  ««  hAg  -  (hAg)B 

‘  ah  h  '  G>  ‘a*  JT — rs - 


(279a) 


According  to  sec.  1,  equations  (8)  and  (9),  this  formula  is  equivalent  to  the 
following  relation  between  the  corresponding  harmonics  of  degree  n: 

('  Th  h!  =  (G^  +  R  (hAg)"  *  (279b) 

so  that  by  (199) 

Ag„*  =  Ag„  +n(|Ag^  +(G,)b  .  (280) 

^  'n 

The  sum  of  the  first  two  terms  on  the  right-hand  side  is  equal  to  the  n-th  degree 
harmonic  of  the  quantity 

Ag  +  n|  Ag  =(\  +  n^Ag  . 

For  lower  degrees  n  (up  to  n  =  5,  say),  the  second  term  is  of  the  order  of 
h/R  and  is  consequently  negligible  as  a  planar  approximation.  For  higher 
degrees  n  this  is  not  true,  because  of  the  large  factor  n.  (Since  a  factor  such 
as  n  always  corresponds  to  some  kind  of  differentiation,  this  is  closely  related 
to  the  fact  that  the  elevation  as  such  may  be  negligible,  but  not  so  its  horizontal 
derivative,  the  inclination.) 

Hence,  for  lower-degree  harmonics, 

Ag*  =  (Ag  +  Gj),  (281) 

or 

("i^f  *0  s  (G^)a  •  (282) 

10" 


This  meins  that  the  Mol odens ky  correction  can  be  applied  to  obtain  in  a  simple 
way  harmonics  of  lower  degree.  For  higher  harmonics  (281)  does  not  hold.  If 
we  would  still  wish  to  apply  Ga  in  this  case,  we  would  have  to  use  the  complete 
expression  (280).  This  expression,  however,  involves  the  spherical-harmonic 
expansion  of  the  product  hAg  and  is,  therefore,  impracticable. 

The  relation  (279b)  was  pointed  out  by  (Molodensky  et.  al,  J962b). 

USE  OF  TERRAIN  CORRECTION.  Let  us  consider  solution  Vllb  of 
sec.  14: 

T  JT  (Ag+G')S(<|>)dff-  £AgAh  ,  (283) 

0 

where 


Ah  =  h  -  h.  , 


(284) 


1\,  being  a  mean  elevation  of  the  region  considered.  In  exactly  the  same  way  as 
we  derived  (225)  from  (222)  and  (223),  we  obtain  from  (283)  and  (223)  the  relation 


G,  =  G'  ♦  |  J-J*- ** -(*■**>,  . 


(285a) 


in  fact,  all  we  have  to  do  is  in  (225)  to  replace  -h  5Ag/Sh  by  G'  and  hAg  by  Ah  Ag/2. 
This  is  equivalent  to  the  relation  between  spherical  harmonics 


n 


(Gi  )n  »  (G')„  -  ^  (Ah  Ag)n  . 


(285b) 


Since  Ah  is  at  most  of  the  order  of  h,  we  may  to  this  equation  apply  the  same 
reasoning  by  which  we  deduced  (281)  from  (280),  finding  that  for  lower  harmonics 
Ag*  -  (Ag  +  G')r  ,  (286) 
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because  then 


This  interesting  relation  between  the  lower  harmonics  (only  those!)  of  three  so 
different  quantities  was  derived  by  Pellinen  (1962,  1964)  in  different  ways. 

The  main  importance  of  the  relation  (287)  rests  in  the  fact  that  G' 
may  in  many  cases  be  identified  with  the  conventional  terrain  correction, 
according  to  (262),  which  can  be  computed  from  the  topography  only,  no 
gravity  anomalies  being  needed.  In  view  of  the  as  yet  imperfect  global  gravity 
coverage  this  furnishes  a  convenient  means  of  estimating  the  correction  (287) 
to  the  lower  harmonics  corresponding  to  the  linear  approximation;  for  a 
practical  estimate  see  sec.  18. 

FIRST-DEGREE  HARMONIC.  From  (278)  follows  that  the  first-degree 
harmonic  Ag*  must  be  zero;  otherwise  the  corresponding  harmonic  Tt  would 
be  infinite.  According  to  (287), 

Ag*  -  (Ag  *  (Ag  +  G^j  +  (Ag  +  G')a  *  O  .  (288) 

This  condition  is  satisfied  by  the  gravity  anomaly  Ag.  It  may  be  shown  that  for 
any  of  the  integral  equation  of  Part  A,  except  (169),  to  be  solvable  a  condition 
equivalent  to  (288)  must  be  fulfilled. 
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17.  THE  EXTERNAL  GRAVITY  FIELD 
If  the  gravity  anomalies  are  given  on  a  level  surface,  then  the  external 
gravity  field  can  be  computed  by  well-known  spherical  formulas  such  as  (200) 
(Heiskanen  and  Moritz,  1967,  chap.  6).  Since  the  free-air  anomalies  Ag  refer 
to  the  physical  surface  of  the  earth,  the  most  natural  way  is  to  reduce  them  to 
a  level  surface  such  as  to  sea  level,  to  obtain  free  air  anomalies  at  sea  level 

Ag*  -  Ag  -  h  (289) 

according  to  sec.  9.  This  expression  for  Ag*  is  valid  to  a  linear  approximation; 
higher-order  approximations  may  be  found  by  an  iterative  solution  of  Bjerhammar's 
integral  equation  by  (174  a,  b). 

Then  the  disturbing  potential  outside  the  earth  is  given  by  (200): 

T(rP,  A,  X)  «  ^  JJ  Ag*  S  (rp,  0,  R)  d a,  (290) 

a 

where 


r„  *  R  +  h„  , 

hp  being  the  elevation  above  sea  level  of  the  point  P  at  which  T  is  to  be  com¬ 
puted.  By  (135)  and  (136)  we  have 

qo  2n  +  1  RP 

S(rp,  4>,  R)  =  l  —r  — n  P„  (cos  4> )  ;  (291) 

n  i  rp 

a  closed  expression  is  obtained  from  (136),  on  replacing  r  by  R. 

The  gravity  disturbance  vector 

—4 

6  =  grad  T  (292) 
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is  found  by  forming  the  partial  derivatives  of  (290)  in  the  usual  way  (Hirvonen  and 
Moritz,  1963;  Heiskanen  and  Moritz,  1967,  chap.  6). 

In  certain  cases,  such  as  for  computation  of  gravity  anomalies  at  flight 
elevations,  it  may  be  of  advantage  to  reduce  the  surface  anomalies  Ag  to  a 
level  surface  at  a  mean  elevation  of  the  area  under  consideration,  rather  than 
reducing  them  to  sea  level.  Reduction  to  point  level  is  theoretically  possible, 
but  less  convenient  for  practical  application. 

Instead  of  analytical  continuation  to  some  level  surface  we  may  also  use 
the  Molodensky  correction  G1 .  In  this  case  Ag*,  eq.  (289),  is  replaced  by 
Ag  +  G1 ,  and  instead  of  (290)  we  have 

T (rp ,  0,  X)  =  ^  P*  (Ag  +  Gt)  S  (rp,  i/i,  r)  da  ,  (293) 

a 

where 

S(rpI  if,,  r)  =  Z  ~~  rr+T  pn  (cos  if,)  (294) 

n=  2  n  1  rp 

in  agreement  with  (135),  equation  (136)  providing  a  closed  expression. 

This  follows  from  Brovar’s  method  expressed  by  equations  (135)  to  (147) 
of  sec.  6*  The  difference  between  (147)  and  (293)  is  only  that  (293)  suppresses 
the  zero-degree  harmonic  of  T,  which  seems  to  have  become  customary  practice 
in  geodesy,  while  (147)  retains  it. 

The  sum  Ag  +  GT  in  (293)  is  by  (137)  only  a  linear  approximation  to  X  sec  /? 
(because  dS  =  R2  do  sec  (9);  higher  approximations  for  X  are  obtained  by  an 
.ierative  solution  of  (142)  by  (143a,  b). 
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Whereas  the  solution  (290)  admits  of  a  simple  geometrical  interpretation, 
being  essentially  a  spherical  formula,  such  as  interpretation  is  not  possible  for 
(293),  since 

r  =  R  +  h  (295) 

is  not  a  constant,  but  a  variable  defending  on  the  topographic  elevation  h.  It  is 
the  spherical  approximation  to  the  radius  vector  of  the  earth's  physical  surface 
(Fig.  4b).  For  practical  application,  the  variability  of  r  constitutes  a  slight 
disadvantage  of  (293)  as  compared  to  (290). 

As  for  other  methods,  the  direct  application  of  Molodensky's  surface 
layer  (sec.  5)  and  the  use  of  Green's  identities,  similar  to  the  way  in  which 
the  integral  equations  of  sec.  2  have  been  obtained,  leads  to  complicated  and 
impractical  formulas  (Brovar,  1963b;  Moritz,  1965b). 

APPLICATION  TO  HEIGHT  ANOMALIES  AND  DEFLECTIONS  OF  THE 


VERTICAL.  The  formulas  just  given  for  the  external  gravity  field  are  applicable 
down  to  the  earth's  surface  and  also  on  the  earth's  surface  itself.  The  height 
anomaly  (  is  then  obtained  by  dividing  the  anomalous  potential  T  by  mean 
gravity  G,  and  the  components  £,  77  of  the  deflection  of  the  vertical  are 
obtained  as  the  horizontal  derivatives  of  f. 

Hence  (290)  and  its  horizontal  derivatives  give  the  formulas  of 
Bjerhammar 


^5  JJA8*s<rp'  *•  "Xfc. 

J 


3S (r„ ,  d>,  R)  r  cos  on 
*4  \  sin*  I**  ’ 


(296a) 
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(296b) 


and  from  (293)  we  obtain  the  formulas 


C  - 


R? 

4trG 


rr  X  sec/3  s  (rp ,  ji,  r)  da  , 
a 


(297a) 


IT  X  secjS 
Or 


a  S  (r„,  r) 

a  ib 


j  COB  <y*i  .  X  COS  Q 
l  sin  ry  J  ^  "  G 


Equations  (296a,  b)  and  (297a,  b)  are  written  in  such  a  way  that  are 
also  valid  for  higher  approximations  to  Ag  and  X,  which  may  be  obtained 
iteratively  using  (174a,  b)  and  (143a,  b),  as  mentioned  above.  The  quantity  r 
is  expressed  in  terms  of  the  elevation  by  (295),  and  since  the  computation  point 
P  lies  now  on  the  earth's  surface,  we  hav 


rp  *  R  +  h,  ,  (298) 

■here  hp  is  the  topographic  elevation  at  P. 

The  solution  (296a,  b)  corresponds  to  free-air  reduction  to  sea  level 
and  was  given  by  Bjerhammar  (1964);  eq.  (296a)  is  identical  with  eq.  (201)  of 
sec.  9. 

The  solution  (297a,  b)  was  published  for  the  zero-order  approximation 
X  cos  /3  =  Ag,  ftt  =<XD  =  /3a  by  Hirvonen  (1960).  Higher  approximations  to  X  were 
found  analytically  by  Brovar  (see  sec.  6),  of  which  the  linear  approximation 

X  sec  0  *  Ag  +  Gj  (299) 

was  used  in  (293). 

The  second  term  on  the  right-hand  side  of  (297b)  is  due  to  the  dis¬ 
continuity  of  the  normal  derivative  of  the  function  (137)  on  the  physical  surface 
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of  the  earth  S.  As  we  have  seen  in  sec.  6,  this  discontinuity  along  the  surface 
normal  has  the  magnitude  -  X.  In  (140)  we  have  used  the  vertical  component 

-  X  cos  0  of  this  discontinuity;  now  the  (negative)  horizontal  components 

-  X  cos  8  tan  0l  and  -  X  cos  0  tan  ft,  figure  in  (297b),  since 


I  9T  1  ST 

*  '  ■  G^r  •  ”  ■  -  o  Sf 

and,  according  to  (Moritz,  1964,  sec.  2.4),  the  unit  normal  vector  n  is  given 
by 


n  =  (  -  cos  0  tan  0X ,  -  cos  8tan  0St  cos  8)  . 

The  functions  S(rp ,  <£,r)  and  ?iS  (rp ,  ft,  r)/dif>  have  been  given  a 

computationally  convenient  form  by  Hirvonen  (1960).  He  sets 

_r_  R+h 
rp  R  +  h,  * 


(300) 


Then  we  have 


Rs  S(  rp ,  tl>,  r)  =  rpt®  [*  ~  +  1  -  3D  -  tcos^(^5  +  3 


in- 


1  +  D  - 1  cos 


3S(rp,  j,  r) 

r)  ill 


,  ,r2  6  _  D-l  +  tcosih 

t?  srn^i  73  +  -  +  3— -  8  - 
L  D3  D  D  sur  ill 

1  +  D-tcos  jfr 


-  3  '’n 


J  l 


2  /j 

(301a) 


(301b) 


where 


D2  =  1  -  2tcos$  + t?  . 

These  expressions  may  also  be  used  in  (296a,  b)  if  t  is  understood  to 
be  R/rp  instead  of  (300). 
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It  is  evident  that  the  system  {296a,  b)  is  considerably  simpler  than 
(297a,  b).  The  quantity  r  is  variable,  whereas  R  is  a  constant;  above  all, 
however,  the  expression  (297b)  for  the  deflection  of  the  vertical  explicitly 
contains  the  inclination  of  the  terrain,  which  is  undesirable  as  we  have  seen 
in  sec.  15.  hi  fact,  the  linear  approximation  of  (297b)  leads  to  (269),  in  the 
same  way  as  (296b)  corresponds  to  (271)  or  (272). 

Since  neither  the  inclination  nor  the  horizontal  derivatives  of  Ag  enter 
in  (296b),  this  expression  shares  the  accuracy  advantage  of  (272);  in  addition, 
(296b)  is  rigorous  in  the  sense  that  it  is  not  restricted  to  the  linear  approxi¬ 
mation.  On  the  other  hand,  the  functions  S{  rp ,  R)  and  %  S(rp ,  ft,  R)/^  j Jj 
are  more  complicated  than  the  simple  Stokes  and  Vening  Meinesz  functions, 
so  that  the  main  practical  significance  of  (296a,  b)  is  its  suitability  for  higher 
approximations. 

The  major  disadvantage  of  the  system  (297a,  b),  corresponding  to  the 
method  of  Hirvonen  and  Brovar,  is  the  inclination  term  in  (297b).  This  term 
would  be  missing  if  the  deflection  components  £  and  r\  were  computed  outside 
the  earth;  the  computation  point  could  be  as  close  to  the  earth's  surface  S  as 
we  like,  as  long  as  it  does  not  coincide  with  S  .  Hence  we  are  tempted  to  try 
instead  of  (297b)  the  formula 


0  “  4,%-  fjx~g  3515^-3  a  a,  , 


where  €  is  an  arbitrarily  small,  but  nonzero,  positive  number;  we  could,  for 


113 


instance,  chooee  c  ■  1  cm.  This  is  certainly  possible,  but  unfortunately  the 
simplification,  as  compared  to  (297b),  is  only  apparent:  if  we  try  to  evaluate 
the  effect  of  the  innermost  zone  in  (302),  we  shall  find  expressions  which 
contain  the  Inclination  after  all,  and  even  in  much  the  same  way  as  (297b)  does. 
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18.  CONCLUSIONS  AND  RECOMMENDATIONS 
in  the  preceding  sections  we  have  tried  to  draw  conclusions  as  to  the 
suitability  of  various  linear  solutions  for  practical  use.  We  may  summarize 
the  results  in  the  following  list. 

Gradient  Molodensky 

Solution  Solution 

height  anomalies  excellent  (I  b)  excellent  (HI  a) 

deflections  of  vertical  excellent  (D  II)  good  (D  III') 
spherical  |  low  excellent  good 

harmonics  high  excellent  poor 

external  field  excellent  good 

The  designations  I  b,  III  a,  DO,  and  D  HI'  refer  to  classification  of 

sections  14  and  15.  The  other  solutions  of  sec.  14  can  be  ruled  out 

practically  because  they  are  less  simple,  except  the  Pellinen  type  solution  VII  b, 

which  will  be  considered  later. 

According  to  the  above  table,  the  best  all-purpose  method  is  the  gradient 
solution.  As  downward  continuation,  or  free-air  reduction,  to  sea  level  it  is 
ideal  for  those  quantities  that  are  computed  on  a  global  scale  such  as  spherical 
harmonics  and  height  anomalies  (or,  which  is  basically  the  same,  geoidal  heights). 
For  deflections  of  the  vertical  a  slight  modification,  reduction  to  point  level,  is 
appropriate.  Here  the  sea-level  anomalies  Ag*  are  not  used  directly,  which 
gives  to  the  computation  of  deflections  of  the  vertical  a  slightly  special  character, 
but  this  is  also  in  conformity  with  the  fact  that  deflections  are  not  usually  computed 
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on  a  global  scale.  Anyhow,  the  vertical  gradient  9Ag/^h  is  used  for  free-air 
reduction  both  to  sea  level  and  to  point  level. 

Computing  this  gradient  in  other  ways  than  (b)  is  not  recommended  for 
general  use. 

Molodensky's  solution  HI,  the  correction  GT  being  computed  by  (a),  is  the 
most  straightforward  expression  for  the  height  anomaly,  although  the  practical 
evaluation  of  Gt  is  more  problematical  than  that  of  dAg/dh  because  GT  is 
strongly  affected  by  the  irregularities  of  topography.  To  get  goexi  results  for 
the  deflections  of  the  vertical,  the  modified  expression  D  ITT  should  be  used  rather 
than  Molodensky's  original  formula  D  HI.  Molodensky's  method  is  of  great  importance, 
but  on  the  whole  the  gradient  solution,  being  simpler  and  more  versatile,  appears 
to  be  preferable  in  practice. 

The  solution  VII  b,  using  free-air  anomalies  combined  with  a  term  G' 
which  is  essentially  the  terrain  correction,  is  fairly  well  suited  for  computing 
height  anomalies,  deflections  of  the  vertical,  and  lower-degree  spherical  harmonics. 

It  is,  however,  net  so  universally  applicable  as  the  other  linear  solutions  since  it 
presupposes  a  strong  correlation  between  free-air  anomalies  and  elevation  (see 
sec.  13).  The  main  importance  of  this  solution  rests  in  the  fact  that  it  furnishes 
a  link  with  conventional  methods.  Since  the  terrain  correction  oan  be  computed 
without  needing  gravity,  estimates  of  G'  are  easily  obtained.  In  addition,  we 
have  seen  in  sections  13  and  15  that  the  uss  of  Ag  +  G#,  free-air  anomalies 
modified  by  adding  the  terrain  effect,  in  Stokes'  and  Vening  Meinesz'  integrals 
without  further  corrections  furnishes  height  anomalies  and  deflections  of  the 
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vertical  at  the  physical  surface  of  the  earth  in  the  sense  of  Molodensky,  at 
least  to  a  better  approximation  than  if  genuine  free-air  anomalies  Ag  were  used. 
This  result  is  quite  gratifying  since  in  the  past  Stokes'  and  Vening  Meinesz' 
integrals  have  sometimes  been  evaluated  using  free-air  anomalies  conventionally 
"corrected"  by  the  terrain  effect;  it  is,  however,  surprising  that  the  results 
are  not  geoidal  heights  and  deflections  of  the  vertical  at  sea  level,  as  one  might 
assume,  but  height  anomalies  and  deflections  at  ground  level. 

This  solution  is  not  suited  for  computing  the  external  gravity  field  and 
spherical  harmonics  of  higher  degree.  Therefore,  there  is  no  reason  for  using 
it  on  a  larger  scale. 

PRACTICAL  SIGNIFICANCE  OF  LINEAR  CORRECTIONS.  The  solutions 
discussed  here  consist  of  spherical  formulas  such  as  Stokes'  integral,  plus 
small  correction  terms.  The  question  arises  as  to  the  practical  significance 
of  these  corrections.  Studies  of  mathematical  models,  such  as  given  in 
(Arnold,  1660),  (Molodenskii  etal.,  1962a)  and  (Molodensky  et  al. ,  1962b), 
are  useful  for  understanding  extreme  situations,  but  they  may  tend  to  give  an 
exaggerated  picture.  For  this  reason  they  should  be  complemented  by  test 
computations  in  selected  areas. 

As  for  height  anomalies,  Arnold  (1960)  has  obtained  a  correction  of  -0.2  m 
for  Mt.  Blftnc,  elevation  4807  m.  This  small  value,  which  is  probably  less  than 
the  error  due  to  the  spherical  approximation  (sec.  3),  seems  to  indicate  that  for 
height  anomalies  the  refinements  of  Stokes'  formula  are  In  general  of  little 
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significance,  although  we  shall  subsequently  see  that  there  is  also  evidence  to 
the  contrary. 

hi  any  case,  matters  are  entirely  different  for  deflections  of  the  vertical, 
ft  should  be  carefully  kept  in  mind  that  0. 2"  in  the  deflection  correspond  to  a 
linear  displacement  of  6  meters,  and  that  the  effort  necessary  to  obtain  the 
deflections  to  an  accuracy  of  0.2”  is  incomparably  greater  than  that  required 
for  the  equivalent  accuracy  of  6  meters  in  the  height  anomaly.  Fig.  11  indicates 
that  small  deformations  (corresponding  to  small  corrections  to  £)  may  give 
rise  to  considerable  changes  of  inclination  (corresponding  to  corrections  to  £ 
and  tj).  Consequently,  it  is  in  the  case  of  deflections  of  the  vertical  that  the 
refinements  of  Molodenskv's  theory  are  of  real  practical  significance.  Except 
in  flat  country,  the  linear  corrections  to  Vening  Meinesz'  formula  are  indis¬ 
pensable  if  an  accuracy  of  +  0. 2",  corresponding  to  an  accuracy  of  +  10  m  for 
the  absolute  position  in  space,  is  aimed  at.  Again,  the  reader  is  referred  to  the 
examples  of  (Arnold,  1960). 

As  for  spherical  harmonics,  especially  those  of  lower  degree,  it  might 
be  expected  that  the  effect  of  the  linear  corrections  is  insignificant,  still  more 
so  than  in  the  case  of  the  height  anomalies.  The  estimates  of  Pellinen  (1962), 
however,  indicate  that  this  effect  is  very  considerable,  of  the  order  of  15  -  20  4- 
of  the  harmonic  coefficients  themselves.  This  surprising  result  seems  to  imply 
that  also  height  anomalies  may  be  noticeably  affected  by  these  corrections.  This 
point  would  probably  deserve  closer  examination. 
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ADEQUACY  OF  THE  LINEAR  APPROXIMATION.  The  opposite  problem 


is  whether  the  linear  approximation  is  practically  sufficient  or  whether  highci 
approximations  should  also  be  considered.  This  question,  just  as  the  preceding 
one,  is  actually  beyond  the  scope  of  the  present  report;  still,  a  few  arguments 
may  be  offered.  It  all  depends  on  whether  the  inclination  may  be  considered 
small  or  not.  An  inclination  ft  of  45°,  so  that  tan  8=1,  is  obviously  not  small, 
but  this  is  certainly  an  exception.  The  average  inclinations  are  usually  small 
even  in  mountainous  areas;  some  kind  of  averaging  is  inherent  in  any  topographic 
map  and  in  any  map  of  gravity  anomalies,  however  dense  the  gravity  survey  may 
be.  Still,  the  fact  that  the  gradient  solution  is  not  directly  affected  by  the  terrain 
inclination  may  be  adduced  as  an  argument  for  preferring  it  to  a  solution  such 
as  Molodensky’s,  in  which  G,  strongly  depends  on  the  irregularities  of  the 
terrain. 

Certainly,  formulas  for  higher  approximations  are  available.  However, 
the  question  as  to  their  practical  applicability  and  to  the  genuineness  of  the 
results  so  obtained  comes  up.  Higher  approximations  correspond  to  the  use  of 
higher  derivatives  of  the  gravity  anomaly  field  and  of  the  topography.  It  is  well 
known  thit  errors  of  a  given  function  are  increasingly  magnified  by  successive 
differentiations.  An  example  well-known  from  geophysical  prospecting  is  the 
second  derivative  ^Ag/^h2 ,  where  different  methods  may  yield  results  that 
diverge  by  more  than  100  *  . 
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These  computational  difficulties  endanger  much  more  than  theoretical 
difficulties  the  use  of  higher  approximations.  Such  theoretical  difficulties  are, 
for  instance,  the  Impossibility  of  rigorous  downward  continuation  to  sea  level 
and  the  corresponding  mathematical  instability  of  Bjerhammar's  integral 
equation  (sec.  7),  or  question  of  the  analytical  convergence  of  the  series 
expressing  the  complete  solution  of  Molodensky's  integral  equation  (Molodenskii 
etal.,  1962,  pp.  122-3),  which  apparently  has  not  yet  been  decided.  But  over 
these  problems,  interesting  as  they  are  from  a  mathematical  point  of  view,  we 
should  not  forget  that  we  are  looking  for  numerical  accuracy,  mathematical  rigor 
being  only  a  tool  to  attain  this  accuracy.  (It  is  well  known  from  mathematics  that 
analytically  divergent,  "asymptotic,  "  series  yield  perfectly  useful  numerical 
results ! ) 

In  keeping  with  this,  close  attention  must  be  paid  to  the  effect  of  inter¬ 
polation  errors,  which  are  present  even  with  a  dense  gravity  net.  In  this  respect, 
conventional  methods  of  gravity  reduction  may  have  advantages:  isostatic  and 
Bouguer  anomalies  are  smoother  and  easier  to  interpolate  than  free-air 
anomalies.  Gravity  reduction  can  indeed  be  incorporated  in  the  "new  theory, " 
but  this  topic  will  not  be  considered  here. 

RECOMMENDATIONS.  Based  on  the  preceding  arguments,  the  gradient 
solution  is  suggested  for  practical  application.  For  computation  of  height 
anomalies,  spherical  harmonics,  and  the  external  gravity  field,  downward 
continuation  to  sea  le/el  is  appropriate,  using  formulas  such  as  (204)  and  (290); 
the  sea-level  anomaly  Ag*  is  computed  by  (199)  and  (177).  For  deflections  of 


the  vertical  the  formula  (272),  which  corresponds  to  reduction  to  point  level,  is 
suggested.  To  combine  these  two  different  free-air  reductions,  to  seg  level  and 
to  point  level,  in  a  manner  suitable  for  large-scale  automatic  processing,  it  is 
suggested  to  compute  and  store; 

surface  free -air  anomaly  Ag, 

correction  -  h  , 
dh 

sea-level  free-air  anomaly  Ag*  , 

5AjI 

anomalous  gradient  — r  ; 

o  h 

the  latter  being  needed  for  the  deflections  of  the  vertical  explicitly,  that  is  not 
only  in  the  form  of  the  correction  -  (dAg/dh)  h. 

If  a  higher  approximation  is  desired,  the  sea-level  anomalies  Ag*  may  be 
computed  by  the  iterative  solution  (174  a,  b);  then  the  height  anomalies  and 
deflections  of  the  vertical  are  obtained  by  (296a,  b).  As  a  matter  of  fapt,  these 
formulas  for  £,  and  r\  can  also  be  applied  if  Ag*  is  computed  by  thf  linear 
approximation  (199). 

These  recommendations  are  based  on  purely  theoretical  considerations 
and  are  consequently  of  a  somewhat  preliminary  nature,  although  they  appear 
plausible.  These  theoretical  considerations  should  therefore  be  supplemented 
by  comparative  computer  studies. 


122 


■■■  rnwwH  fpj 


APPENDIX  I 


PLANAR  APPROXIMATIONS 


PLANAR  APPROXIMATION  TO  STOKES'  FUNCTION.  As  a  planar 
approximation,  the  first  term  on  the  right-hand  side  of  (16)  was  neglected 
to  get  (176).  In  the  same  way  we  may  neglect  the  corresponding  terms  in 
equations  (10)  and  (11),  whereby  they  become  identical.  To  this  extent, 


also  their  inverses  (1)  and  (5)  are  equivalent,  so  that  the  planar  approximation 
to  Stokes'  function  is 


S(iM  = 


SIMPLIFICATION  OF  (103).  In  sec.  4  we  have  found  the  solution  (103): 


T  ■  ^  jr  (ag + p> + 1® ") s  <*>  do  * • 


We  shall  now  show  that  the  term  3q,  /2R  between  the  parentheses  is  negligible  as 


a  planar  approximation. 

According  to  (97)  we  may  set 


-qx  *  A  +  B, 


(305a) 


%  =  A  +  2B  , 


(305b) 


where 


A  -£  <h-MTda 

a 


(306a) 


(306b) 


Let  us  first  consider  the  term  A.  Its  effect  on  T  is 

8T‘  '  £  IT(£><*»*-  <3°7> 

a 

If  we  compare  A  with  GT ,  oq.  (226),  we  see  that  the  only  difference  is  that  Ag 

in  Gj  is  replaced  by  T  in  A.  By  comparing  (222)  and  (223)  we  find 

£  .rjO,S(*,d»  .  hS«,)da-hAg; 

a  a 

whence,  on  replacing  Ag  by  T, 

£  ;j  A  s(4 >  da  -  -  £  jr  fh  h  S  w  <te  -  h  T , 
a  a 

so  that  (307)  becomes 

fT*  *  £  II  a  O * ' f)  s<*>  * ‘ft  T  •  <308> 

0 

Now  the  integrand  is  of  the  order  of  (h/R)  Ag  which  is  negligible  with  respect 
to  the  gravity  anomaly  Ag,  which  is  the  main  term  of  the  integrand  of  (304), 
and  the  last  term  of  (308)  is  of  the  order  of  (h/R)  T,  which  is  negligible  with 
repect  to  T.  Hence  the  complete  effect  6TA  is  negligible  as  a  planar  approxi¬ 
mation. 

As  for  B,  the  comparison  of  (306a)  and  (306b)  shows  that  A  and  B  are 
of  the  same  magnitude  since  sin  $  dh/d^i  is  of  the  same  order  of  magnitude  as 
h  -  hp  for  small  .  Therefore  we  may  conclude  that  the  effect  fiTB  will  be 
of  the  same  order  of  magnitude  as  6Tt  and  will  consequently  be  also  negligible 
as  a  planar  approximation. 

The  weak  spot  in  this  argument  in  that  the  above-mentioned  equality  of 
orders  of  magnitude  of  sin  $  Bh/rty  and  h  -  h,  holds  for  small  ;0  (of  the  order 
of  a  few  kilometers)  only.  It  is  therefore  desirable  to  have  an  independent  check, 
although  the  computations  will  be  somewhat  lengthy. 
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4  t 


We  consider  the  effect  of  qt ,  which  by  (305a)  is  the  combined  effect  of  A  and  B: 
ATC  =  6Ta  +  6Tb  ;  (309) 

we  shall  set 


■  £ J7  rKh  *  ^ -sin*^')Td®' 

a 


(310) 


By  (220)  this  is  equivalent  to 


„  R2  r*  r»  sin  l/)  „  ,  .  *  T  . 

0  =  *  *> 
a 


(310*) 


or 


C  =  D  +  E 
where,  by  (38), 

D  -  -  h-  27  ff  -rf  *  ■ h  £  !!  5  (T  •  r  )  • 


(311) 


R2  ,.p  sin  ib  .  dT  . 

~  U  „  3  “  • 


2  V 


(312a) 


(312b) 


The  transformation  of  (312a)  by  means  of  (33)  yields 

o 

Since  AgT  =  *Ag/^h  by  (180),  this  becomes 
D  -  -  h  f-  ff  — 

2?r  JJ  fe  \  *h  J 

o 

and  since  this  integral  is  the  planar  inverse  of  (176),  the  expression  for  D 


reduces  to 


D  -  h  Ag  ; 


(313) 
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hence  the  effect  on  T  is  in  analogy  to  (307)  given  by 


»T°  -  i;  If  M**s<*>* 


As  for  (312b}  the  effect  of  E  is 

"« -  i  xr  H  s<*>  * 

o 


i-|  jx  ■■!#)-». 


Now  (312b)  corresponds  to  (231),  h  replacing  Ag,  and  T  replacing  h.  Hence 


the  integral 


4ir  JJ  ES(ri)da, 


which  corresponds  to  (232),  is  by  (240)  equal  to 

*  iTh  r57  s<*>  *  ■ 

a 

so  that  the  effect  of  £  becomes 


f rp  R  pn  3  h  .  , 

STt  Tn  !}  2R 

o 


Adding  (314)  and  (315)  according  to  (311)  we  find  the  effect  of  C  =  -qt  to 


_R  pc  3h  /  A  A  *  T 
4 it  U  2R  ‘g  R>^ 


Ag  +  r>j6  ^  S  ^  ^  * 


Since  all  components  of  the  "gravity  disturbance  vector"  grad  T  will  be  of  the 

same  order  of  magnitude,  the  horizontal  component 

ftT 
R *ib 


will  be  of  the  order  of  the  vertical  component 


a  r 


Ag  + 


2T 

R 


and  consequently  of  Ag  itself.  Hence  the  effect  of  C,  eq.  (316),  is  negligible 
as  a  planar  approximation  for  the  same  reason  as  (308)  is. 

As  a  by-product  we  find  the  effect  of  B  by  subtracting  (308)  from  (316) 


to  be 


, —  K  p  p  3h  /'  -  T  2  T  N  ,  3h 

B  =  4ir  ?!  2R  C  R*i h  R  -  (,j)  60  +  2R  T  ’ 


2T 


3h 


(317) 


whose  planar  approximation  is  likewise  zero. 
Hence  (304)  reduces  in  fact  to 


T  =  4ff  TJ  (Ag  +  P.JStjiJda  +  q,  , 
a 


(318) 


which  is  identical  witi  (104). 
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APPENDIX  n 


MEAN  CURVATURE 

We  shall  now  derive  the  formula  (191)  for  the  difference  of  mean  curvature 
between  two  surfaces  S  and  S ,  whose  distance,  measured  along  the  normal  to  S, 
is  £. 

A  formula  for  a  similar  purpose  is  eq.  (8)  of  (Weatherbum,  1930,  p.  173); 
however,  it  is  laborious  to  adapt  it  to  our  present  case.  Therefore  we  shall 
proceed  in  a  more  direct  manner. 

Denote  the  position  vector  of  S  by  x  and  its  normal  vector  by  n  .  Then 
the  relation 

AgX  =  -  2Jn  (319) 

holds,  where  A2  is  the  surface  Laplacian  defined  in  sec.  1.  This  formula  is 
equivalent  to  expressions  given  in  (Weatherbum,  1927,  p.  231)  and  (McConnell, 
1931,  p.  203,  example  3),  only  the  notations  for  mean  curvature  aie  different. 

Our  J  corresponds  to  the  negative  of  McConnell's  H  and  of  Weatherbum's  J/2. 

We  shall  use  a  tensor  notation  equivalent  to  McConnell's,  except  that  we 
are  using  ordinary  instead  of  Greek  indices  for  surface  tensors,  and  a  notation 
such  as  x  ,  without  any  indices,  for  space  vectors.  Then  (319)  may  be  written 

2Jn  =  -  a'  3  ^  }  =  -  a1  >  (x, ,  -  {  ,k.  }  x* "),  (320) 
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j 


_AJ 


where 


a, ,  =  x,  •  x, 


(321a) 


••s  the  first  fundamental  tensor  of  the  surface,  and  x,  and  x,  ,  denote  partial 
derivatives  with  respect  to  the  surface  coordinates.  The  second  and  third 


fundamental  tensors  are,  as  usually, 


i  "  ”  *  nj  *  x,  j  *  n 


c, ,  =  a  ’  n,  =  -  n,  ,  •  n  , 


(321b) 

(321c) 


with  the  relation 


c,  ,  =  ak  1  b,  k  b, ,  ( 

connecting  the  three  fundamental  tensors.  If  K  is  the  Gaussian  curvature  of 
the  suriace,  then 

a!  is,,  =  b‘  '  b,  j  *  4Ja  -  2K  .  ( 

These  relations  may  be  found  in  any  textbook  on  tensor  analysis,  such  as 
(McConnell,  1931). 

— * 

The  position  vector  x  of  S  is  evidently  given  by 

x  «  x  +  £  n  .  ( 

By  (324)  and  (321  a)  we  have  for  the  first  fundamental  tensor  of  S 

a«.j  "  a!<  '  2Cbi  j  •  I 

retaining  only  terms  linear  in  £ .  The  contravariant  fundamental  tent-or  is 
given  by 

a1 3  =  a1  J  +  2£b*3  ,  I 
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since  by  (325)  and  (326) 

a,  3  a3k  *  a,3  aJk  *  ft,k  . 

We  further  set 

Q- 

where  we  need  not  evaluate  the  linear  residuals  r!3  . 

Multiplying  (320)  by  n  and  taking  (321b)  into  account  we  find 

2  J  -  -  a<  3  b,  j  .  (328) 

Let  us  now  multiply  the  equation  corresponding  to  (320), 

2 J n  =  -  a  ' 3  x.  , 

by  n  .  Then,  since  as  a  linear  approximation 
n  •  n  =  cos  0=1, 

we  have 

— ♦ 

2J  =  -  a* 3  xt  >3  •  n  .  (329) 

By  differentiating  (324)  we  find 


~  - 

#  lr  \  — • 

X!,j  =  Xl, 

-  { ij}  x‘ 

"  *1  3  + 

Cl,  n  +  Ci  "j  +C,n,  +Cn 

{?,) 

■  K*C*n  +  CM  -  A,k 

so  that 

X1  ,  ,  *  11  =  kjj  +  Cl  ,  "  C  C1  ,  ”{ij}  C* 
=  +  Ct ,  j  “  Cci  j  • 

130 


By  means  of  this  equation  and  of  (326)  we  obtain  from  (329) 


2J  =  -  (a1  '  +  2£b!  ')  (b,  .  +  f\  ,  ,  “Cc,  .) 

=  2  J  -  2  £  b!  5  b, ,  -  a*  ■  £,  f .  +Ca‘-  c,  ,  , 

so  that 

J  =  J  -  (2  J2  -  K)  C  -  |  A-  C  • 

Thus  the  derivation  of  (191)  is  finished. 

This  equation,  without  derivation,  was  given  in  (Moritz,  1962). 
A  similar  formula  was  published  by  Marussi  (1957). 
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